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Abstract 

We study the Maxwell-Bloch model, which describes the propagation of a laser through a 
material and the associated interaction between laser and matter (polarization of the atoms 
pg ■ through light propagation, photon emission and absorption, etc.). The laser field is described 

through Maxwell's equations, a classical equation, while matter is represented at a quantum 
level and satisfies a quantum Liouville equation known as the Bloch model. Coupling between 
laser and matter is described through a quadratic source term in both equations. The model 
also takes into account partial relaxation effects, namely the trend of matter to return to its 
natural thermodynamic equilibrium. The whole system involves 6+N (N + l)/2 unknowns, the 
six-dimensional electromagnetic field plus the N (N + l)/2 unknowns describing the state of 
matter, where N is the number of atomic energy levels of the considered material. 

We consider at once a high-frequency and weak coupling situation, in the general case of 
anisotropic electromagnetic fields that are subject to diffraction. Degenerate energy levels are 
allowed. The whole system is stiff and involves strong nonlinearities. 

We show the convergence to a nonstiff, nonlinear, coupled Schrodinger-Boltzmann model, 
involving 3+iV unknowns. The electromagnetic field is eventually described through its enve- 
lope, one unknown vector in C 3 . It satisfies a Schrddinger equation that takes into account 
propagation and diffraction of light inside the material. Matter on the other hand is described 
through a ^-dimensional vector describing the occupation numbers of each atomic level. It 
satisfies a Boltzmann equation that describes the jumps of the electrons between the various 
atomic energy levels, as induced by the interaction with light. The rate of exchange between 
the atomic levels is proportional to the intensity of the laser field. The whole system is the 
physically natural nonlinear model. 

In order to provide an important and explicit example, we completely analyze the specific 
(two dimensional) Transverse Magnetic case, for which formulae turn out to be simpler. 

Technically speaking, our analysis does not enter the usual mathematical framework of 
geometric optics: it is more singular, and requires an ad hoc Ansatz. 
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1 Introduction 

Maxwell-Bloch systems are of common use in Laser Physics (see the textbooks [8], [TO], [IS], [54] , 
[56] . [3T]). They modelize the evolution of an electromagnetic field, treated classically, and coupled 
with an ensemble of identical atoms, which in turn are described by a quantum density matrix. This 
model is relevant when atoms are far from the ionization energy (to possess discrete energy levels), 
while they have sufficiently low density and the laser field is strong enough (which allows to describe 
the field classically while matter is described in a quantum way - see |16j). 



In the Maxwell-Bloch model, the electromagnetic field satisfies Maxwell's equations, whose un- 
knowns are the electric and magnetic fields E G M. 3 and B G R 3 , where E = E(t,x,y,z) and 
B = B(t,x,y,z) and t G M is time while (x,y,z) G M 3 are the space coordinates. Matter is de- 
scribed through a Bloch equation, whose unknown is the density matrix p = p(t, x, y,z), a quantum 
variable which describes the atomic state at (t,x,y,z). We consider that the atoms only visit the 
iV lowest energy levels. The latter are the N first eigenstates of the free material system, in the 
absence of fields. In this basis, the density matrix p = p(t, x, y, z) is an iV x iV matrix, for each value 
of (t,x,y,z). The diagonal entries p(t,x,y)(n,n) (called populations) give the proportion of matter 
that lies in the n-th energy level (n — 1, . . . , N), while the off-diagonal entries p(t, x, y, z)(n,p) with 
n ^ p, (called coherences) give the correlation between levels n and p. The complete Maxwell-Bloch 
system takes into account the coupling between the laser field and the atoms via terms that are 
quadratic, proportional to p x E, and which describe polarization of matter due to laser propagation. 

We study the high frequency and weak coupling behaviour of the Maxwell-Bloch system, a situa- 
tion in which the typical frequencies of the field and of the atoms' oscillations are large and possibly 
resonate, while the strength of laser-matter coupling is small. 

For large frequencies, the electromagnetic field is expected to be asymptotically solution to a 
nonlinear Schrodinger equation. This is the paraxial approximation. We refer to [TS] and [2S], as 
well as [30] on these matters, when the sole laser field propagates (no coupling with matter). When 
matter is actually coupled to the field, we refer to [9j. Here a high frequency Maxwell-Bloch system 
is studied both physically and mathematically, for atoms that only possess 3 non-degenerate energy 
levels. The analysis leads to a Schrodinger-Bloch approximation of the original system, in a spirit 
similar to the present paper. 

The weak coupling behaviour of matter is a bit more delicate to handle: to have a clean limit, 
one needs to take thermodynamic fluctuations into account. For this reason we introduce, in a 
standard fashion, phenomenological relaxation operators in the original Bloch system. These impose 
a rapid decay of coherences, as well as a quick return to equilibrium of populations. We refer to 
[3] for mathematical properties of the relaxation operators that are natural in this context. Due to 
the relaxation effects, it is expected that Bloch's equation is asymptotic to a Boltzmann equation, 
sometimes called "Einstein's rate equation" (c/. |33j . [1]). It describes how the atoms jump between 
the various energy levels under the action of the external field. When the driving high frequency 
field is given (and the Bloch model is thus linear), we refer to the papers [6] and [7], which study 
the actual convergence of Bloch's equation to a Boltzmann model in the weak coupling regime. In 
that case, a formula is found for the transition rates involved in the limiting Einstein rate equation, 
which coincides with the one formally obtained in the Physics literature. Note however that the 
question studied in [6] and [7] is a linear problem, and proofs strongly use ODE averaging techniques 
as well as the positivity of relaxation operators (features that the present text does not share). We 
also mention [5], where similar asymptotics are treated both for quantum and classical models. We 
stress finally that many other works deal with the rigorous derivation of Boltzmann like equations 
from (usually linear) models describing the interaction of waves/particles with external media. A 
non- convergence result is given in [H] and |15j . Convergence in the case of an electron in a periodic 
box is studied in (TTJ, [12], [13], while the case of an electron in a random medium is addressed in 
[22] . [29] , [35] . [35] . [10] - see also [35] for a semi-classical approach. In a nonlinear context, a partial 
result is obtained in [2J. 

The above formal discussion suggests, in the present case, that the high-frequency Maxwell system 
goes to a Schrodinger model for the envelope of the field, while the weakly coupled Bloch system 
supposedly goes to a Boltzmann equation describing the jumps of electrons between the atomic levels. 



This is the program we rigorously develop in the present paper. We fully prove convergence of 
the coupled Maxwell-Bloch system to a coupled Schrodinger-Boltzmann model. We also prove that 
the rate of exchange between the energy levels is proportional to the laser's intensity. In doing so 
we recover the physically relevant model. Our approach mainly uses three-scales geometric optics, 
yet in a more singular context where the partial relaxation effects impose a specific treatment of 
coherences. 

2 Presentation of the results 

2.1 The model 

The Maxwell-Bloch system, whose unknowns are the electric field E = E(t, x, y, z) G M 3 , the magnetic 
field B = B(t, x, y, z) G M 3 , together with the density matrix p = p(t, x, y, z) G C JVxAf (the space of 
N x N complex matrices), reads 

d t B + cur\E = 0, (1) 

d t E - curl B = -d t P, with P = Tr (I», (2) 

dtp = -i[n-E-r, P ] + Q(j>), (3) 

where Q(p) = W (jp d - 7p od . (4) 

In these equations, curl is the usual curl operator on vector fields in IR 3 , T as well as Q are given 
matrices in C NxN , 7 > is a given positive constant, and whenever A and B are matrices in C NxN , 
the trace Tr (A) denotes the usual trace of A while the bracket [A, B] denotes the usual commutator 
between matrices 

[A, B}= AB- BA. 

The term Q(p) = Q(p)(t, x, y, z) is the so-called relaxation matrix , an N x N matrix for each value of 
(t,x,y,z). Its definition involves W G C NxN , a given matrix with nonnegative entries W(n,m) > 0, 
while pa and p 0( j denote the diagonal respectively off-diagonal parts of the density matrix p (they 
correspond respectively to the populations and the coherences). They are N x N matrices defined, 
for each value of (t, x, y, z), by their entries 

p d (t,x,y, z) (n, p) = p(t,x,y,z) (n, p)l[n=p], 
pod(t,x,y,z)(n,p) = p(t,x,y,z)(n,p)l[n^p\. 

Equation (J3J) also uses the following notation, valid thoughout the present text : given any matrix 
A with nonnegative entries we set 



N 

A (j p d (n, n) = ^ l A (k, n) p d {k, k) - A(n, k) p d {n, n)] 

fe=i 
A ft Pd( n , p) = when n ^ p. 



(5) 



The meaning of operator Q(p) in ([3]) is the following. The term — 7p d induces exponential relaxation 
to zero for the coherences, while the term W fj p d acts on the populations only, and induces exponential 
relaxation of the populations towards some thermodynamical equilibrium that depends on the values 
of the VF(n,p)'s. As in conventional kinetic theory of gases, relation ([5]) asserts that along time 
evolution, atoms may leave with probability W(k, n) the fcth eigenstate to populate the nth eigenstate 



(this is the so-called gain term Y2k=i W(k, n) pa{k, k)), while some atoms may conversely leave with 
probability W(n, k) the nth state to populate some other /cth state (this is the so-called loss term 
-T,k=i w ( n , k )pd(n,n)). 

Theoretically, the Maxwell-Bloch system needs to be supplemented with the Ampere and Faraday 
laws, 

divS = 0, div(£ + P) = 0. (6) 

These constraints ([6]) are anyhow transported as soon as they are satisfied by the initial data, hence 
we shall skip them in the sequel. 

Note that the above equations are readily given in the convenient dimensionless form that suits 
our purpose. The precise scaling under study is discussed later. 

We now comment on these equations, and on all involved quantities. 

The density matrix p(t, x, y, z) is Hermitian and positive. It describes the state of matter at 
(t,x,y,z). 

The constant matrix Q is the free Hamiltonian of the material system, written in the natural 
eigenbasis. It is a fixed physical constant associated with the considered atomic species. It reads 

n = diag(w(l),...,w(JV)), (7) 

where < w(l) < • • • < u(N) are the atomic energies. For later convenience, we readily introduce 
the differences between energy levels, as 

u)(n,k) = u)(n) — u)(k). (8) 

The constant matrix T is called the dipolar operator. It is a hermitian matrix. It has the value 
-^atomic x 7atomic, where iV atomic is the number of atoms per unit volume, while 7 a tomic is a fixed physical 
constant (a matrix) associated with the considered atomic species. 

The entries of so-called dipolar momentum E ■ V are defined for any m,n = 1, . . . , N, as (E ■ 
T){t,x,y, z)(m,n) = E(t,x,y, z) ■ T(m,n), where a ■ b denotes the componentwise product of two 
vectors in C 3 , namely 

(E ■ T)(t, x, y, z)(m, n) = E x (t, x, y, z)T(m, n) x + E y (t, x, y, z)T(m, n) y + E z (t, x, y, z)T(m, n) z , 

and the subscripts x (resp. y, resp. z), denote the x (resp. y, resp. z) components of the relevant 
vectors (note the absence of complex conjugation). 

The non-negative transition coefficients W(n, k) > are known as the Pauli coefficients, see |10j . 
They satisfy a micro-reversibility relation at temperature T, i.e. W(n, k) = W(k, n) exp (u(k, n)/T). 
In that perspective the relaxation operator Q(p) translates the fact that the atoms tend to relax 
towards the thermodynamical equilibrium given by p d = and Pd{n, n) = exp (—u(n)/T). Both the 
off-diagonal relaxation term 7 and the Pauli coefficients W(n, k) are physical data. Contrary to T 
or Q, their relation with given physical constants attached with the specific atomic species at hand 
is unclear. 

The Bloch equation ([3]) relies on the so-called dipolar approximation : the coupling between light 
and matter is taken into account through the simplest [E-T, p] term, which is quadratic (proportional 
to p x E) and local (it only depends on the value of E and p and the same point (£, x, y, z)). This 
approximation implicitely assumes that the wavelength of the field is larger than the typical spatial 
extension of the atom, hence can be taken constant over the whole domain occupied by each given 
atom. 



2.2 The scaling 

Let us introduce the physical scales in the model, which transform all constants (T, Q, 7, W) and 
unknowns (E, B, p) into quantities of order one in the regime we wish to study. 

Firstly, concerning the time variations, we study a high frequency regime. Calling 1/e the di- 
mensionless parameter measuring typical values of the frequencies, i.e. the ratio between the time 
scale of observations and the time scale of the variations of E, B and p, all time derivatives d t then 
become edt in the scaled version of (pQ)- @j. 

Secondly, we want to study a situation where constructive interference occurs between the time 
oscillations of the electromagnetic field and the ones of the atom (giving rise to atomic absorption 
and emission of photons). For that reason, the time variations of E, B, and p should all take place 
at similar frequencies. For that reason Q is naturally a quantity of order one in the scaled model. 

Thirdly, the typical strength of the coupling between light and matter is entirely determined 
by the physical constant I\ Since weak coupling is realized when the polarization operator has an 
effect of the order one on the chosen time scale, it turns out that we need to prescribe T = 0(y/~e). 
Mathematically, this means we shall replace T by \J~eY in the scaled model. Indeed, the electric dipole 
momentum then becomes \E-(yfeY)\ 2 = 0(e) which, integrated over macroscopic time scales 0(1/ e), 
results in an energy contribution of order one as desired. Physically, since T = N atornic x 7 at0 mic where 
-^atomic is the number of atoms per unit volume and 7 a tomic is the polarizability of one atom, this 
means that the atomic density is here tuned to be of order 0(^/e), so as to observe order one effect 
of coupling over this scale. 

Fourthly, concerning the space variables, the hyperbolic nature of the Maxwell equations ([I])-© 
suggests to rescale space so that d x , d y , d z become ed x , ed y , and ed z , respectively, due to the scaling 
d t *->■ ed t and to finite propagation speed: the laser visits space scales 0(1/ e) over time scales 0(1/ e). 
This simple scaling would in fact be easily described along lines similar to the present analysis. We 
wish to investigate physically richer situations where diffraction occurs. To this end, as in [9], [19] 
and [20], we impose anisotropy and introduce a third scale: on the one hand, we choose one direction 
of propagation, say the x direction, meaning that at the microscopic scale, the fields vary with x 
but not with (y,z); on the other hand, we restrict our attention to fields that slowly vary over the 
scale 1/e in the x direction (or, in other words, that have spatial extension 1/e in x), while they 
slowly vary over the scale 1/y/e in the (y, z) direction (or, in other words, that have spatial extension 
1/y/s in (y, z)). The typical shape of the laser beam is thus that of a "light cigar", as in [17]. In 
macroscopic scales, this provides fields that strongly oscillate at frequency 1/e in x (and only in this 
direction), while they have support of size 1 in x, and support of size \fe in (y, z). All this imposes 
to rescale d x as ed x , and (d y ,d z ) as (yfed y) \fed z ) in the above equations. We stress that this is a 
definite choice of shape of the kind of laser beams we wish to study. It is also a choice of polarization: 
these beams are shot in the x direction only. 

Fifthly, consider relaxations Q(p). In order for the diagonal relaxation W$p& to have an 0(1) 
effect at times t of order 0(1), we take the coefficients in W of size e, and write eW(n, k) instead of 
W(n, k). The off-diagonal relaxation — 7p d on the other hand is supposed to have a much shorter 
time scale, which we choose to be of order 0(e). This corresponds to an off-diagonal relaxation 
— 7 Pod that remains unsealed. In fact, there is no theoretical description of the relaxation time as 
a function of e. In |6] and [7J, the off-diagonal relaxation is scaled as — 7£ M p d with an extra free 
parameter p that is constrained to satisfy < p < 1 (to have a clean limit, off-diagonal relaxation 
should be strong enough with respect to the chosen time scale). For technical reasons, in these 
papers, p is actually restricted to < p < 1/4. In the present paper, the formal analysis could be 
performed for < p < 1, yet our main stability result (Theorem 15.21) requires the stronger constraint 
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/i = 0. 



After rescaling all variables and physical constants accordingly, the system (jT|)- (J4j) becomes 
( d t B £ + curl £ E £ = 0, 

d t E £ - curl e B £ = -J= Tr (r fi 7 C £ ) - *Tr (r[£ e • T, C £ + N e ]) -^Tr (TW$N £ ) } 



< d t C £ = --Q^C £ + -^\E £ -T,C £ + N £ 
e y/e 

d t N £ = ^ r [E £ ■ T, C £ ] d + W ft N e . 



Od; 



(9) 



Here and in the sequel, we adopt for convenience the notation C £ and N e for coherences p Q d and 
populations pd, respectively. Besides, throughout the sequel, subscripts "d" and "od" shall always 
refer to diagonal and off-diagonal parts of the considered N x N matrices. Lastly, we also denote by 
curl £ the curl operator associated with our scaling, i. e. 

curl e E = ( —i=d y E z - —d z E y , —=d z E x - d x E z , d x E y - —=d y E x J , (10) 

and fi 7 is a shorthand notation for ad^ — ij, 

fi 7 C = [Q, C] - i7C, i.e. (n y C)(n,p)=uj(n,p)C(n,p)-ijC{n,p). (11) 



2.3 Description of the results 

Main result: profiles, separation of scales, and obtention of an approximate solution. 

Maxwell-Bloch's system (Q is a nonlinear hyperbolic symmetric system, singular in e, that we may 
write symbolically 



L (d t , d x , -±=d y , -^d z \ W = F £ (V £ ). 



Its unknown is 



U £ = (u £ ,C £ ,N 



£ ,C £ ,N £ ), 



where u e stands for the sole electromagnetic field 

u e = {B £ ,E £ ). 

We are interested in solving a Cauchy problem associated with Q for initial data that are smooth, 
but high frequency, of order 0(1/ s). The typical difficulty is to ensure existence of the whole family 
(U £ ) ee ]o, Eo ] on some time interval [0, £*], independent of e. 

This enters the framework of 3-scales diffractive optics (see [IS] , [30], and the surveys [27], |2T]). 
To this end, our analysis uses a WKB analysis, based on profiles that we expand in successive powers 
of e. An important original point yet is, we consider amplitudes that are even higher (namely 0(1)) 
than the ones allowed by transparency properties (namely 0(e) - see |28j). We show these larger 
amplitudes are eventually compensated by the partial relaxation: remember the sole equation on 
the coherences C £ carries the relaxation term — r )C £ /e in (Q. We refer to remark l3~2l below on this 
important point. 



We start with a given a wavevector 

k= (k\...,k d ) eR d , 
and initial data of the form 

Uf ni (x, y,z) = TP (x } y, z, *fj + 6* L y, z, ^ . (12) 

The profiles U^(x, y, -2,6*0), £ e (^, y, z, 9q) 6 C°°(IR x T d ) that are periodic in 6* 6 M d and we assume 
that 

in every Sobolev space H". We also assume that U^ is submitted to polarization conditions, namely 
that the electromagnetic field u° has the usual directional constraints, see (155]) . and that some 
preferred components C°(m,n) vanish, see (136]) . 

We build for all e > an approximate solution U app to (J2J), defined on some time interval [0, £*], 
and with initial values U°(x. y, 2, kx/e). Our construction uses a leading profile U° and two correctors 
U 1 , U 2 , and U app is of the form 



2 
U a PP (*, x, y,z) = Y^ v^ J U J (t, x, y, z, T, a, O , #1 



(13) 

T=t/ v / e i cr=7t/e,(e ,6»i)=(A;x/e,-fct/£) 



Correctors U 1 , U 2 , are introduced so as to ensure consistancy of this Ansatz, namely 

L [d u d x , -Ld y , ±d)j U app - F £ (U app ) = 0(y/i). (14) 

One main point in our analysis is the separation of scales, and the crucial introduction of the fast 
variables T, a and Q\ in ( fl3j) . Explanations on this point are postponed to the next paragraph. 

Once the approximate solution is constructed, we prove Theorem 15.21 It asserts that, for e > 
small enough, there is a unique solution U e to system (j§j) with initial value Uf ni , which is well 
approximated by U app , namely 

V/i e N 3 , \\d^ z (U £ - U app ) || Loo([0A]xR3) ->0. (15) 

On top of that, we are able to completely describe the dynamics of the dominant term U°. It provides 
an o(l) approximation of the original dynamics of U e . The function U° satisfies a coupled, nonlinear, 
Schrodinger-Bloch system, which we describe later. 

This result is achieved via a singular system method (see j25j), where the unknown is a profile 
(thus with non-singular initial data). The difficulty comes from the "supersingular" nature of the 
system, and the fact that relaxations are only partial ones (acting on a part of the dependent variables 
only). This problem is overcome thanks to the structure of the approximate solution. 

Describing the asymptotic dynamics (1) - oscillations and initial layer. In the chosen 
Ansatz ( 115]) . we assume that each profile LP may be decomposed into modes U^ K , as 

W(t, x, y, z, T,a,9)= J2 J2 W <? ex P (*' a ' 9 ) ex P ("^ " ( 16 ) 
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This defines the quantities U^ ,K = XJ^ K (t,x,y,z,T), for which a £ 1? d is seen as a Fourier mode, 
while k £ N is seen as an exponentially decaying mode. We introduce in passing and for later 
convenience the notation 



U J a = ^Urexp(- K a), (17) 



kGN 



a function of (t, x, y, z, T, a), as well as 



T> K = J2 U « K ex P (»«^)> ( 18 ) 

agZ 2d 

a function of (t,x,y, z,T,9). For k > 0, the decomposition (fT6l) encodes an exponential mode in 
the a = jt/e variable, representative of the initial layer induced by off-diagonal relaxations, while 
Fourier modes a = (ao,ai) reflect oscillations in the variable (0 ,9i) = (kx/e, —kt/e). 

This choice of Ansatz is motivated by the following. The oscillations in kx/e are anyhow present 
in the initial data, and the hyperbolic feature of the equations ensures they are propagated into 
oscillations in the variable (kx/e,— kt/e). Nonlinear interaction of waves next makes sure that all 
harmonics are created along time evolution in (Q. To be more precise, we show that all oscillations 
that are characteristic for the Maxwell-Bloch system are propagated. They are given, for fields and 
populations, by the characteristic set 

C := C + U C_ U C , 

where we define 

C± = {a = (a , «i) £ Z 2d \{0} | «! = ±« } and C = {a = {a , ot x ) £ Z M \{0} | ot x = 0} , (19) 

corresponding respectively to the waves propagating to the left and to the right, and to purely 
spatial oscillations. The reader should be cautious about the fact that the Fourier mode a = is not 
considered as a part of this characteristic set. These sets essentially correspond to the characteristic 
variety of the linear Maxwell part of the equations. For coherences, we show the characteristic 
frequencies are those which resonate with some transition energy u(m,n). They are given, for any 
wavenumber k £ Z d , by the resonant set 

U{k) = {{m, n, a , on) £ {1, . . . , iV} 2 x Z 2d | k ■ on = u{m, n)}. (20) 

Describing the asymptotic dynamics (2) - Rectification effects. The rectification phe- 
nomenon is the creation of non-oscillating terms through the nonlinear interaction of oscillating 
terms. This phenomenon is typical of quadratic systems such as (JHJ). 

It is the reason for our introduction of the intermediate time scale T = t/y/e in ( TT6l) . This 
scale captures the evolution of the system between the macroscopic length \x\ = 0(1) and the 
wavelength 0(e). The first need for the intermediate scale T is, the non-oscillating terms induce a 
secular growth of the necessary corrector terms (see [31J). In that circumstance, it turns out that 
imposing a sublinearity condition, 

-||U j || L oo — ► for j = l,2, (21) 

ensures smallness of correctors in (fl3l) . namely y/eXJ^t, x, y, z, t/y/e, o~, 9) = o(l) when e — > as 
desired (j = 1,2). A second reason for the introduction of the scale T is that, already in the less 
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singular regime of diffractive optics, the leading profile U° obeys some linear, constant coefficients 
hyperbolic evolution system with respect to time T, and the T-sublinearity condition for correctors 
provides a unique way of determining U° (see [26], [30], and [19], [20] for the variable coefficients case), 
via a detailed analysis of wave interactions at the scale T (see Section H~3"j) . A similar phenomenon 
occurs in the present situation as well. 

More precisely, the above mentioned system at scale T here takes the following form. Separating 
the average part and the oscillating part as (we use the notation introduced in (fl~6]) - ([T8]) ) 

u° = Uq + u° sc , (which defines u° sc = u° — u°), 

we prove that u(j satisfies a system of the form 

d T u° + M 2 (0,d y ,d z )u° = 0, (22) 

where M2(0, d y , d z ) is a maira-coefficient differential operator of size 3x3. Eventually, our analysis 
shows that (122]) has to be solved together with the one giving the corrector uj, namely 

d T ul + M 2 (0, d y , d z )ul = -M x {d t , 8 X ) u°, (23) 

where Mi(dt,d x ) is another matrix differential operator. We establish that equation (123]) possesses 
a unique solution uj provided T-sublinearity of the right-hand-side is imposed. It turns out that the 
same "secular growth" analysis is also necessary for u° sc and N°. 

Remark 2.1. In the particular case of profiles that do not depend on T, the coupled system ( 122]) . 
( 123]) is actually overdetermined (see Remark \6.1\) . Solvability is only recovered provided the matrix 
operators Mi(d t , d x ) and M 2 (0, d y , d z ) commute. This very particular and important situation occurs 
in the Transverse Magnetic case (see below), where these operators are all scalar. In the general case, 
rectification enforces the introduction of time T to make the set of profile equations solvable. 

Describing the asymptotic dynamics (3) - the coherences. For coherences, we establish the 
relation 



c o_ c o,i e - CT= I J2 C°/e ia A e-° 
\aez 2d / 



In other words, at dominant order coherences decay as exp(— jt/s). Besides, the coefficients C^' 1 are 
shown to satisfy (see (14"6"]) ) 



V(m, n, a) e K(k), d T C^ n>a = i [£ u ■ T, C^] m ^ a . (24) 

The set lZ{k) is defined in ( 12T)]) . The other components C^ n a , for which (m, n, a) <£ IZ(k), are shown 
to vanish (and E° is determined independently: see the paragraph "Field dynamics" below). 

Now the resolution of ([24]) cannot rely on the above mentioned sublinearity condition, due to 
the fact that a generic solution to (124]) grows like exp(KT) or so at least, for some constant K > 0. 
Similarly, the analogous equations on the correctors C 1 and C 2 cannot be solved using the sublinearity 
condition neither. As a consequence, the evolution of C° remains a priori undetermined. The key 
point now comes from the off-diagonal relaxations. Indeed, we do prove that C 0,1 satisfies a bound 
of the form 
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for some K\, K 2 > 0, so that the product (C ' 1 exp(— a)) \ T _ t/ r cr __ tj has size Oi exp(K 2 t/y/e) x 

x exp(— jt/e) j, hence is negligible when e — > 0, as desired. We notice a posteriori that the lack of 

knowledge in the evolution of C° is harmless, since coherences live during an initial layer of size 0(e) 
only. We also show that the correctors C 1 and C 2 are negligible thanks to estimates in the same 
vein (see Section |4. 3. 2p . 

Describing the asymptotic dynamics (4) - the electromagnetic field. First, we show that 
u° = u 0,0 , so that the electromagnetic field does not undergo the same decay as coherences. Next, 
we prove that u°, conveniently decomposed into its average and its oscillatory part, may be written 

u° = u° i0 + u° >+ + u° _ + u t ° ime + < ace , (25) 

* v ' * „ ' 

— u° =u° 

where we define the purely spatial oscillations, and temporal oscillationqj, as 

o = V u° e ia ' e u°. = V u° e ia - e 

"space / j u o ° i "time / , u a ° 

aSCo a&C+UC- 

The three mean terms Uq , Uq + and u[j _ are defined in the course of the analysis. The above 
decomposition entails the fact that we have u° = whenever a (£ Cq U C + U C_ . This is a polarization 
condition, as usual in optics. Besides, there are algebraic projectors 7To, 7r + and 7r_, such that for any 
a G C t , we have 7r t u° = u° (t = 0, +, 

i. 

space i 



We also show that functions Uq , u° ce , u° + , Uq_ do not depend on the slow time t, i.e. 



«9 t < = 0, u s ° pace = 0, ftug i± = 0, (26) 

and they satisfy, from ([22]) . 

d T u 0fi = 0, 9 T u° P ace = 0, (dr ± ^l + Dl) ug, ± = 0. (27) 

The latter is a two-dimensional (non local) transport equation at time scale T. The time oscillations 

u time satisfy 

\/a&C±, d T u° a = 0, and (28) 

(d t ± d x ) u° a T ^\X = 7r ± (0,^Tr (r(fi 7 (fi 7 -A; • a,)' 1 - l)[E? imc ■ r,N°] Q ) ), (29) 

a nonlinear Schrodinger equation with respect to time t, describing diffraction in the transverse 
variable (y, z), and quadratic coupling between light and matter. Technically, a standard Diophantine 
hypothesis ensures that k ■ a is "not to small" in (l2"9~j) - see Hypothesis 13.31 

Describing the asymptotic dynamics (5) - the populations. As for the electromagnetic field, 
the dominant term N° = N°'° in the populations does not decay in time at first order. In addition, 
we show that N° only possesses spatial oscillations, thanks to the identity 

N° = K + N° pacc . 



Ho be accurate, these oscillations, namely u° ime , involve both time and space variables - we neertheless keep the 
denomination "time" for this part of the oscillations. 
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This is a polarization property. Furthermore, N° does not depend on the intermediate time T, 

«9 T N° = 0. (30) 

Lastly, at the macroscopic time scale t, populations evolve according to a Boltzmann-type equation, 
with transition rates that are the sum of the linear rates W (see ©), and of nonlinear rates that 
depend quadratically on E° and change with the frequency 9. Precisely, we get the relation 



aN° = w a n° 



(e° osc + e°) ■ r, (i n, -k ■ fcj^pk + e°) ■ r, n 



01 
od 



d, space 



(E° osc + E° 0>0 ) • T, (« Q, -k ■ de^-'KE ^ + E° 0>0 ) ■ T, N°] od . (31) 

J d,0 



We refer to (JHH]) below. In practice, this equation is of the form 

a'+a"=a 

for some Pauli coefficients W(o/,a",E^,,,E^„_g) that depend on the Fourier indices a' G Z M , a" G 
Z M , (3 G Z 2d , as well as linearly on each variable E®, and E Q a „_^. 

An illustrative example : the Transverse Magnetic (TM) case. We detail here the above 
analysis in the classical TM case (see [M]), for which formulae become more explicit, and the geometry 
of the problem is a bit simpler. The TM case is a 2-dimensional case, for which in addition fields 
have the particular polarization 



B = B(t,x,y)= \B y \, E = E(t,x,y) 




Firstly, we give the explicit profile equations in this specific situation. Secondly, we show that the 
introduction of the intermediate time T is not necessary in that case (see Remark 12. II) . Finally, we 
perform the analysis for prepared data, i. e. when coherences vanish at leading order. Eventually, we 
show that the approximation is better than in (|15[) . in that Theorem 16.31 asserts 

Vfi G N 2 , H^, (U £ - U^ pp ) IlL-dp^xB?) = 0{y/e). 

2.4 Outline of the paper 

In Section [3], we describe the Ansatz (IT31 . the choice of scales and phases. We also discuss the needed 
Diophantine hypothesis 13.31 on the wave numbers (k l , . . . , k d ) G M. d involved in ( 1T2"|) . In Section HI 
we proceed with the construction of an approximate solution U| pp = XJ e (t,x,y, z,T,9) that is 
consistant with Maxwell-Bloch's system ([9]). Sections 14.11 and 14.21 contain the fast scale analysis, 
i.e. the analysis in the variable 9. This is a Fourier analysis, leading to the usual characteristic sets 
and group velocity of geometric optics. Section 14.31 is devoted to the intermediate scale analysis of 
the resulting profile equations at the time scale T. The sublinearity condition (12 ip is involved (in 
Section f4.3.ip . that allows to treat wave interactions. We also perform (in Section |4.3.2[) the ad hoc 
analysis for the coherences C| pp , which have exponential growth in T and violate the sublinearity 
condition. In Section 14.41 we eventually solve the Cauchy problem for profile equations in all variables 
t, T, and 9. To complete the analysis, Section [5] provides the proof of our main result, namely that 
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U| pp indeed approximates the true solution U e to the Maxwell-Bloch system (jSj), see equation ffTBT) 
above. The precise form of our result is given in Theorem 15.21 Lastly, Section [6] goes again through 
the whole analysis, yet in the simpler Transverse Magnetic case, and for prepared initial data. The 
corresponding result is given in Theorem 16.31 

Our main theorems are Theorem 15.21 and Theorem 16.31 

3 Formulating the Ansatz 

We solve the Cauchy problem associated with ([9]) (with unknown U e = (u £ , C e , N e ), u £ = (B £ , E 6 )) 
for rapidly oscillating initial data of the form 

Uf ni O, y, z) = Y^ 0K X > ?/' z ) ex P (* G 9 ' k ) -) ' ( 32 ) 

/3ez d 

where the wave vector k is k = (k 1 , . . . , k d ) G M. d for some d G N*, its coordinates are chosen Q- 
independent, and the scalar product (3-k in (1321) denotes (3-k = (3 1 k 1 + - ■ ■ + (3 d k d . The Q-independence 
of fc's coordinates ensures that for (3 G Z d , relation (3 ■ k = holds if and only if (3 = 0. The vector k 
collects the independent oscillations (in x) carried by the initial signal (B e , E £ ,p £ ) ini (x,y, z), while 
the integer /3 J G Z numbers the various harmonics corresponding to the phase W x/e (J = 1, . . . , d). 
We choose an initial signal that carries all possible harmonics. This harmless choice is motivated 
by the fact that the original, nonlinear, Maxwell-Bloch system anyhow generates nonlinear wave 
interaction, which implies that any initial oscillation creates the associated harmonics along wave 
propagation. 

The off-diagonal relaxation term — jC 6 /e in (Q enforces exponential decay of coherences C £ , so 
that, at first order, we expect only the fields u £ = (B e , E e ) and the populations N e to propagate. 
Hence space and time oscillations are expected to be generated from the initial spatial ones only 
through Maxwell's equations and through the equation for the populations. In a similar spirit, 
since our scaling postulates that variations in the transverse variables (y, z) only occur at the scale 
y 7 ^ (and not e), we also expect that oscillations will not occur in the (y,z) variables. In other 
words, the relevant Maxwell equations for oscillations are expected to be ID in the x direction. The 
corresponding characteristic variety, in Fourier variables, reads 

C M axwell U Copulations = { (t, £) G 1R 2 \{0} | T = Or T 2 = £ 2 } . 

Eventually, propagation of space and time oscillations are expected to occur via the collection of 
phases 

k? x — io° t, with to 3 = or uj 3 ' — ± k 3 (J — 1, . . . , d). 

More precisely, the relevant oscillations in our analysis are expected to be all the harmonics (3 ■ (k x — 
cut), as the multi-index /3 runs in Z d , and the j-th coordinate of u is either 0, k 3 or — k 3 . An instant 
of reflexion shows that this ensemble coincides with the collection of phases (a\ ■ k) x — («o • k)t as 
oq and a.\ run in Z d . This motivates the following 

Notation 3.1. For any multi-index a = (a 1 , . . . , a 2d ) G Z M ; we decompose 

a := (a , ati), with a = (a 1 , . . . , a d ) G Z d , and ai = (a d+1 , . . . , a 2d ) G Z d . 
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Accordingly, any (smooth enough) function V(0) defined over the torus T 2d may be decomposed in 
Fourier series as 

v(0) = J2 Ya ex p (* a • e "> = J2 v ° ex p (» ( a ° • ^o + «i • ^i)) • 

aGZ 2d a£Z 2d 

With this notation, resonances are identified by the characteristic set C = Co U C + U C_ , where 
Co, C + , and C_ have been defined in (fl9j) . 



Due to the fast decay of coherences, another set of phases needs to be introduced. Namely, the 
off-diagonal relaxation term — jC 6 /e in §§§ leads to introduce the complex phases exp(— jt/e) as well 
as all its harmonics exp(— Kjt/e), where k£N. 

With all these considerations in mind, we can now define the approximate solution we seek, as 
Uapp = « pp , C| pp , N| pp ), < pp = (Bl pp , El pp ), where 

(33) 



T=t/^e,a=-Yt/e,(8o,di)=(kx,-kt)/£ 



Ul pp (t,x,y,z) = ^2y/e 3 \J J (t,x,y, z,T,a,9 Q ,9i 

3=0 

and for any j = 0, 1, 2, and any a > 0, T > 0, 9 E T 2d , we set 

U 3 (t,x,y,z,T,(T,9) = ^ J2U£"(*,a;,y, z,T) exp (ia ■ 9) exp(-Kcr). 

oez 2d «eN 

Remark 3.2. (discussion of the chosen profiles). 

(%) Coming back to the original scales of the problem, we may represent the solution U e under the 
form 

U £ (t, x, y, z) = U £ (t, x, y/ey, y/ez), 

where the function U £ has variations in y = y/ey and z = y/ez at scale y/e, and solves with the 
operator curl e replaced with curl, while the associated initial datum is of the form 

lJ^ (x,y,z) = V^ (x,^=,^=)=:U^Jx $ * X 



e y/e/ '-" \ \/e \/e e 



This is the 3-scales setting of [20] (where only the case of quasilinear, non- dispersive systems is 
adressed). It leads to a profile representation of the form 

U e (t, x, y, z) = U £ ( t, x, y, z, — , - 



where <f> = (f>(t,x,y,z) is the collection of the two phases kx and —kt, while ip = ^{t-,x : y,z) is 
a collection of "intermediate phases" (in ( |33l) . U e does not depend explicitly on y,z, so that if) = 
(t,y,z)). When non-oscillating terms are present at first order (XJq / 0), it is in general necessary 
to put in this collection if) an intermediate time t in order to solve the profile equations via the analysis 
of Section \4-3\ This is the reason why the variable T = t] y/e is present in (|33|) (it may be unnecessary: 
see f20\f . Remark 1.4, and the Transverse Magnetic case below, Section® for which if) = y). This 
intermediate time T captures the evolution of intermediate scales y/e, between "macroscopic length" 
(9(1) and wavelength e. 

(ii) The consistancy of the chosen Ansatz requires that y^U 1 <C U° whenever a = jt/s, T = t/ y/e, 
9 = (kx/e, —kt/e). This requirement enforces T-sublinearity of the non- exponentially decaying part 
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(corresponding to k — 0) of the correctors U J ' K (j > I), see (16TJ) . The exponentially decaying 
correctors (k > 1) may be of the same order as the first profile U° A : they anyhow lead to o(y/e) error 
terms, see Section \4.3.2 . 



(Hi) The analogy with 1201/ shows that we could treat the same problem with curved (i.e. nonlinear) 
phases <fi. This situation arises in particular when dealing with inhomogeneous media, where variable 
coefficients (i.e. variable electric and magnetic permittivity) are involved in the Maxwell-Bloch system, 
(iv) We emphasize the fact that in the present work, we deal with large population variations, of order 
one. 

In I2BJ (long time diffraction, p. 248), with no relaxation terms, transparency allows to trans- 
form the original Maxwell-Bloch system, via a change of dependent variables, into a dispersive and 
quadratic system of the form L(ed)U = e 2 f{U), with U = (U^U 11 ), where (B,E,C) = ell 1 , and 
N = e 2 U n . The system is considered over times of order 1/e. Taking relaxations into account, we 
consider here a similar system, at the same space and time scales (oscillations have frequency 1/e, 
and propagation is considered over times of order 1/e), yet in a situation where (B, E, N, C) are 
larger, of order 1, and the quadratic coupling term f(U) is stronger, of size e, namely 

ed t C £ = -i[Sl, C% d + ie[E £ ■ T, C £ ] od + %e[E £ ■ T, N% d - 7 C £ , 
ed t N £ = ie[E £ ■ T, C £ ] d + s 2 W (1 N £ . 

Of course, the stronger quadratic interaction term is balanced by the off-diagonal relaxations: note 
however that relaxation only affects part of the unknowns (the coherences), and one key aspect of our 
analysis precisely relies in the analysis of the stronger interaction term. 

Such large population variations are also considered in /£j|/, without relaxation, yet for waves that 
are polarized in a specific way, and in the case of a three-level Bloch system only (i.e. N = 3 in our 
notation). In that case, additional conservation properties are at hand. On top of that, /2J/ considers 
a weaker coupling term of size e 2 . The limiting system obtained in f9] is of Schrodinger-Bloch type. 

Since small divisor estimates naturally enter the analysis below, we readily formulate the usual 
Diophantine assumption we shall need on the wave-vector k. It will be used in order to invert the 
differential operators acting in the 6 variable on the various (smooth) profiles LP ,re (t, x, y, z, T, 9). We 
refer for instance to [21]. 



Hypothesis 3.3. The wave vector k G 1R S from fl32|) is Diophantine, namely 

3C,a>0, V/3GZ d \{0}, |/3 • k\ > C\f3\- a . 

Remark 3.4. The above assumption is harmless. Indeed, the following fact is well-known. Pick any 
exponent a > d — 1. Then, the set 

{k e R d I 3C > 0, V/3GZ d \{0}, |/3 ■ k\ > C\f3\- a } 

has full measure in M. d . In other words, almost any k G M. d (for the Lebesgue measure) has the 
Diophantine property. 

4 Formal expansions and approximate solution 

Notation 4.1. Denote by M(dt,d x ,d y ,d z ) the order 1 differential operator in Maxwell's equations, 



curl 
— curl 



M(d t , d x , d y , d z )=(_ " ^ )=d t + A x d x + A y d y + A z d z 
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where the Aj 's are 6x6 real symmetric matrices. Set also 

M 1 (d t , d x ) = M(d t , d x , 0, 0) = d t + A x d x , M 2 (d T , d y , d z ) = M(d T , 0, d y , d z ) = d T + A y d y + A z d z . 

With this notation at hand, plugging the Ansatz (133]) into the Maxwell-Bloch system (Q, we get 
that U| pp satisfies the original Maxwell-Bloch equations up to a residual r e that is given as follows 
(here we use the notation (TTT1) for fly) 

Proposition 4.2. The residual 



r £ := L U, d x , -Ld y , ±=d)j Vl pp - F £ (U: pp ) 



has the profile representation 



r £ (t, x, y, z) = n £ (t, x, y, z, T, a, 



T=t/^/e,a='yt/s,(eo,e 1 )=(kx,-kt)/e 



where 



K £ =J2J2 V^V>(*, *, V> ^ T , e ) exp(-«:a) 

j=-2 k>0 



and the first terms (—2 < j < 0) are (see Section [3T71 for the others) 

( Mi(-jK-k-d 01 ,k-de o )u°> K \ 



-2.K 



V 



(i f2 7 — 7K — k ■ dg ± )C 
(-7« - k ■ d ei )N°' K 



0.K 



/ 



/ Afi(- 7 « -k-d 9l ,k- d do ) u 1 ^ + M 2 (<9 T , fy, &) u°- K - (0, i Tr (r fly C°' K )) \ 



-1,K 



(i f2 7 -jK-k- doAC 1 '" + d T C°' K - % [E° ■ T, C° + N 



01 K 

Jod 



V 



-7/c 



fc • dflJN 1 '* + <9 r N°' K - z [E° ■ T, C°]3 



/ 



„0,K 



N 



0lK x 



/ Ml (- 7 «: -k-d 6l ,k- d 9o ) u 2 ' K + M 2 (d T , d y , 8 Z ) u 1 '- \ 

+M l (d t , d x ) u°' K + (o, -i Tr (r f2 7 C 1 '*) + % Tr (r[£° • T, C 

(i fly -jK-k- d 8l )C 2 ' K + d T C 1 ^ + «9 t C°> K 

-i [E° • r, C 1 + N 1 ]^ - % [E l • r, C° + N°]* d 

I (-7« - fc ■ fll )N 2 > K + o^N 1 '" + <9 t N°' K - z [E° ■ T, C 1 ^ - i [E 1 ■ T, C% - W (j N°- K 

With this computation at hand, we wish to construct the first profile U°' K = (u°' K , C°' K , N°' K ), 
in C°([0,t*]t,C°([0,+cx)[T,-H' oo (M 3 x T 2d ))) for some U > 0, together with correctors (u-»,C J ',N-») 
(j = 1,2), in such a way that the residual r e is small (cf. Section [5.2p . More precisely, we shall 
impose r~ 2 = r _1 = r° = 0, and show that this procedure completely determines U°. This is 
obtained by decomposing successively (with increasing j) the equations obtained for the LP's, and 
separating characteristic and noncharacteristic modes (a,n). The necessary linear algebra tools are 
developed in the next paragraph. 
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4.1 Rapid modes and algebraic projectors 

• Tools needed to deal with the electromagnetic field. The Fourier series representation, in 
the 9 variable, of the first equation stemming from [r _2,K = 0] is (according to notation 14. ip 

Mi(-7« - ik ■ a u ik ■ a ) u° a K = 0, V«6 Z 2d . 

For each a, this is a system of linear equations in C 6 , with matrix 

M\(— 7« — ik ■ ai,ik ■ a ) = ik ■ aj. + i(k ■ a )A x , 

which may be singular only if k = 0. To deal with this equation, we need the following classical 
definitions and lemmas 



Definition 4.3. For each a £ Z 2d , let ir a be the orthogonal projection in R 6 onto the kernel of 
Mi(—ik ■ a\,ik ■ ocq),. Denote by M\(—ik ■ a\,ik ■ ao) _1 the inverse of M\(—ik ■ a\,ik ■ ao), when 
restricted to the space orthogonal to its kernel, namely M\[—ik ■ a%,ik ■ ao) -1 acts on Ran(l — ir a ). 
Finally, define the projector II on the space of Fourier series, 



/ I /y Hi qj O • 



n Y, »° e "-° : = E 

\aeZ 2d J a£Z 2d 

A straightforward computation establishes the 

Lemma 4.4. For each a £ Z 2d , the projector TT a is a homogeneous function of a of degree zero. It 
takes a constant (matrix) value tc 1 ' on each component C L of the characteristic set (i = +, — , 0), the 
value one (or identity) for a = 0, and vanishes else. In particular, we have the identity 

n ( Yl u ^\ ■= u ° + Yl K +U * eia ' 6 + J2 *~ u c*J a ' 6 + J2 K°u a e ia - e . (34) 

/ aeC + aeC- aeC 



In any circumstance, for any Fourier series u = 2_. u a e ta ' , we have 

a&Z 2d 

Mi(—k ■ d 6l , k ■ dg ) u = ^=^> Hu — u. 

Next, the following two lemmas are classical for geometric and diffractive optics with smooth 
characteristic varieties ([32], [21], [15] , and also [H] for an elegant unified version). They express 
that the operators acting on the oscillating part of the fields are in diagonal form. 

Lemma 4.5 (group velocity). 

For any Fourier series u = 2_. u a e% ■> containing no mean term, we have, with v(o) = i for 

a£Z 2d \{0} 

a £ C L (l — 0, +, —), the relations 

TlA y Ilu = TlA z Tlu = } and TlAJlu = ^ v(a)n a u a e ia - 9 =: v(D )Tlu. 

a£C 
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Lemma 4.6 (diffraction). 

For any Fourier series u = N u a e ia ' e , containing no mean term, we have 

a£Z 2d \{0} 



UM 2 (0,d y ,d z )M 1 (-k-d ei ,k-d eo )- 1 M 2 (0,d y ,d z )Uu 

^ &(a)A y)Z 7i- a u a e ta ' e =: i&(D e )A ytZ Ilu, 



- i 



isc+uC- 



where a(a) = ±— for a G C± and a(a) = for a G C . 

2k ■ a 

Proof. By the definition of the operators involved, for each a, one has 

TT a M 2 (0, d y , d z ) M x (-k -ai,k- a ) _1 M 2 (0, d y , d z ) ir a = ^ KaAjM^-h ■ a u k ■ a^ 1 A k -n a djd k . 

j,k£{y,z} 

Now, 7i a and M\{—k ■ «i, fc ■ a ) _1 are the evaluations, at £ = (k ■ a , 0, 0), of the spectral projector 
P{t (£),£) resp. of the pseudo-inverse M _1 (t(£), £) of the complete Maxwell symbol M(r(£), £). Here 
and in the sequel, we use the notation 

7"±(0 = =F|ei for a G C t ,± := C L n {A; • ai ^ 0}, t = ±, 

together with 

r(0 =0 for a G C . 

The quantities t±(£) and r(£) are naturally the eigenvalues of the Maxwell symbol M(t(£),£). On 
the other hand, we have the identity, valid for any j, k G {x, y, z}, see [21] or [4"T] . 



<^% 



= P(r(0, A,- M- 1 (r(0, ^ ^(r(0, + ^(0, ^ M" 1 ^), ^ P(r(0, 0- 
Using now the obvious relation 



x ^ (0 = ^ id -WrTu: rDiag (0,1,1) for £= (A; -00,0,0), 

inspection of the five cases a G C Ly ± (t = ±), and a G Co, leads to the lemma. □ 

Lastly, in order to distinguish between propagated and non-propagated parts of the profiles, 
we introduce the following splitting, refering to the "oscillating" part of the fields, to the "time" 
oscillations of the density matrix, or to its "space" oscillations. 

Definition 4.7. For any Fourier series U = ^2 aeZ 2d U a e ia ' e , we set 

TT ■— \~*TT p ia ' e TT — V^ TT p ia ' e TT ■ ■— V^ TT p ia ' e 

w osc • / , w a° ? "-^spacc • / ^a 1 - j "-^timc • / ^n 

a^0 aeCo aeC+UC- 
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• Tools needed to deal with the coherences. For each mode (a, k), the second equation from 

j r -2,K _ gj reac [ s 

Vm, n, (z(w(m, n) - fc • «i) + 7(1 - «))<%»,« = 0. 
Here, only k — 1 is of interest. This justifies the introduction of the resonant set lZ(k) we defined in 

Remark 4.8. In [9], the wave vector k is precisely chosen so that TZ(k) be nonempty. 

• Tools needed to deal with the populations. The populations are scalar variables, and the 
characteristic frequencies are simply the (a, «)'s belonging to C x {0}. No additional tool is needed. 

4.2 Profile equations, fast scale analysis 
4.2.1 The residual r 2 

According to the notation above, the equation [r~ 2 = 0] is equivalent to the polarization conditions 

u° = u°'°, Vk > 0, u°- K = 0, and u° = Ilu , (35) 

C° = C ' 1 e- CT , and V(m, n, a) <£ K{k), C^ Q = 0, (36) 

N° = N - , and N° = N° + N° pace . (37) 

In order not to overweight notation, we shall from now on systematically refer to u° and N° in 
the sequel, keeping in mind they do coincide with u 0,0 and N 0, °. 



4.2.2 The residual r 



-1 



Equations on the field. We first deal with the mode k = 0. Separating the mean term (a = 0) 
and the other Fourier modes, we find 

d T u° + M 2 (0,d y ,d z )u° = 0, (38) 

d T u° 0SC = 0. (39) 

Obtaining the second equation requires to use (135]) and Lemma |4~5| which imply I1M2(0, d y , d z )U = 
and u° sc = IIUq SC . The next order profile u 1 ' is then seen to satisfy 

(1 - nK;° c = -M 1 (-k -d dl ,k- cy- 1 M 2 (0, d y , d z ) u° osc . (40) 

Secondly, when k > 1, the operator M±(— 7K — k ■ de t , k • de ) is invertible, and we get the two 
values 

u 1 ' 1 = Mi(-7K - k ■ d 0l ,k ■ de.Y 1 (0,zTr (r ^C ' 1 )), (41) 

\/k > 1, u 1 ^ = 0. (42) 
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Equations on the coherences. We first note that whenever k^I the operator (i f2 7 — r jK 1 — k-de Q ) 
is invertible. Thanks to (135]) . (|36|) and (|37|) . this gives the two values 

C 1 ' = i (i fi 7 -fc • d ei )- l [E° ■ T, N°] od , (43) 

Vft > 1, C 1 '" = 0. (44) 

When «; = 1, we need to distinguish between resonant and non-resonant triples (m,n,a), to 
obtain 

V(m, n, a) <£ TZ(k), C^ Q = -(u(m, n) - k ■ a^E ■ T, C^ l ] m ^ a , (45) 

V(m, n, a) G fc(*), d T C° m ^ a = i[E°- T, C ' 1 ]^,*- (46) 

Equations on the populations. Here we need to distinguish between the values ft = and ft 7^ 0, 
as well as between the modes a G Cq U {0} and the other Fourier modes. In that way we obtain, 
taking ft = and restricting to modes a G C U {0}, the relation 

<9 T N° = 0, (47) 

together with the following polarization conditions and spectral properties for N 1 (here we use (13710 

Va i C U {0}, N^'° = 0, N 1 ' 1 = i(- 1 -k-de 1 )- 1 [E -T,C°' 1 } d , Vft > 1, N 1 '* = 0. (48) 

4.2.3 The residual r° 

Equations on the field. When ft = 0, we have (using (J131) and (1361) ) 

= Afi(-fc • d 6l , k ■ d 6o ) u 2 - + M 2 (Qr, 9„, <9 Z ) u 1 ' + M x {d u d x ) u°'° (49) 

+ (0, i Tr (T[E°'° • r, N ' ]) - z Tr (Tz ft 7 (z ft 7 -k ■ deJ'^E ' ■ T, N°'°] od )) . 

Again, we separate oscillating and nonoscillating parts in the above equation. In the nonoscillating 
case, the two nonlinear contributions are seen to compensate each other and we obtain 

M 2 (d T , d y , 8 Z ) ul'° = -M 1 (d t , d x ) u°, (50) 

On the other hand, the Il-polarized part of oscillating modes turns out to satisfy the relation, using 
(USD, (BOD, (USD and Lemma SSI and writing M 2 (d T , d y , d z ) = d T + M 2 (0, d y , d z ), 

d T Uulf c = -UM 1 {d t , d x ) u° osc + UM 2 (0, d y , d z ) M 1 {-k -d ei} k- d^)' 1 M 2 (0, d y , d z ) u° sc 

+ n(o,zTr (r^^-i^)- 1 -!)^!,^)). 

Thanks to Lemma [4.51 and Lemma [4. 6[ this is rewritten, using u° sc = ITu^, 
d T nul£ =-(dt + v(D,))u° sc + i &(D e )A y>z u° 08c 

+ n(o,zTr (Vii^i^-k-de^-l^-V^X^y (51) 

Lastly, applying (1 — II) to the oscillating part of equation (j4"9"|) . we also recover (using the fact that 
8t ((1 — II)u ^) = 0, deduced from equation (j4"0l together with the identity <9tu° sc = 0), 

(1 - n)u 2 f c = - Mii-k -d 9l ,k- de Q y l (M 2 (0, d y , d z ) u o ; c + M 1 (d t , d x ) u° sc 

+ (0,z Tr (r(z fi 7 (z ft 7 -k ■ de,)- 1 - l)[E° ■ r,N°] osc )) ). (52) 
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When k — 1, using (T4"Tj) . (1551) . (157j) . we get in the same way 

u 2 > 1 = -M 1 (- 1 -k-d ei: k-de )- 1 [M 2 {d Tj d y7 d z )M 1 (-k-d eiJ k-do )- 1 (0,zTr (T(i ^C - 1 ))) 

- (0, z Tr (r fi 7 C 1 ' 1 )) + i Tr ([£° • I\ C ' 1 ]) 1 . (53) 

Finally, for greater values of k, using (jl"Z|) . ( 155]) . (pHj) . ( 155]) . (j57|) . we have 

Vk > 1, u 2 - K = 0. (54) 

Equations on the coherences. When k / 1, we may invert directly 

c 2 ' K = -(ia r -iK-k- do,)- 1 (d T c l ' K + d t c°' K -i[E°- r, c 1 + n 1 ]^ - % [e 1 ■ r, c° + N°]* d ) , (55) 



which vanishes for « > 2, thanks to (|44j), (1361), (T48|), (157j). 

For k = 1, using (1361) . (1571) . (jlij) . (jlSj) and (j4"2"|) . nonresonant triples (m,n,a) lead to the similar 
formula, 



N 1 

J-' Jm,n,Q 



V(ro,n,aO£fc(AO, C^, a = - (iw(m,n) + 7 - f* ■ air^C^ - i [E° ■ r,^' 1 

(56) 

- i [e 1 ' • r, c - 1 ]™ - 1 [E 1 * 1 • r, N ] m ,n,, x 



whereas for resonant triples we obtain 

V(m, n, a) e K(k), d T C l m ) n ^ a = - d t C° m ] n , a + z[E - 1\ C 1 ' 1 + N 1 ' 1 



+ i [E^ • T, C - 1 ]^ + i [E 1 ' 1 ■ T, N°] m , n>Q . (57) 

Equations on the population. For k = 0, using (1551) . (156|) . (1571) . (1441) and (|48|) . we have 

Va G Co U {0}, 9 T N^° = -d t N° a + t[E°- T, C 1 ' V + W J| N° , (58) 

Va £ C U {0}, N 2 ' = (i fc • ax)" 1 (d^ + «9 4 N° - i [E° ■ T, C 1 ' ]^ - W fj N°) . (59) 

Similarly, k > leads to a polarization relation (using (157)) ). namely 

\/k>0, N 2 ' K = (-fK + k-d 9l )- 1 (dT'N 1 ' K -i[E -T,C 1 }^-i[E 1 -T,CX\ (60) 

which vanishes as soon as k > 2. 

4.3 Profile equations, intermediate scale analysis 

The next step in the analysis consists in obtaining a closed system determining the first profiles 
u°, N° and C°. To achieve this, and in order to ensure consistancy of the Ansatz, we need impose 
that the corrector terms u 1 , N 1 and C 1 are small compared with the first profiles u°, N° and C°. 
Concerning the field and the populations, this means we need impose T-sublinearity of the correctors 
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u 1 ' and N 1 ' (as in [23], [2E], [SO]) while solving equations (150]) . ( 15T]) . ( 155]) . and we shall prescribe 
the following requirement 

VaGN 3 + M i sup R^^N 1 ' )^ — > 0. (61) 

J te[0,t*] T^+oc 

Concerning coherences, i.e. while solving (157]) . we cannot impose the same constraint on C 1 ' 1 . 
However, and as explained before, this corrector produces in the approximate solution a term 
y / £C 1 ' 1 (t, x,y, z, T, 6)\ T= e~ 7t//£ , so that the possible growth in T of C 1,1 is eventually compen- 
sated by the factor exp(— jt/e). This is proved in section H.3.21 

4.3.1 Analysis of fields and populations: average operators 

The key observation for solving equations (150]) . ( )5Tl) . and (158]) on u ' , IIu^, and N ' , while keeping 
T-sublinear solutions u ' , ITu^, and N ' , is that the source terms in these equations have a precise 
stucture in terms of propagation at the intermediate scale T. 

Let us make our point precise. Our analysis is in three steps. 

Firstly, equations (150]) . (I5T|) . and ( 155]) are of the form 

<9 T nu^ c = • • • , a T Nj'° = • • • , d T ul>° + m 2 (o, a„ cy uj'° = • • • , 

where the right-hand sides only depend on the lower order terms u° and N°. Besides, for (77, () G 
M 2 \{0}, the symmetric matrix M 2 (0,/7,C) has the spectral decomposition 

M 2 (0,7/,C)= Yl X k(v,C)Pk(v,0, (62) 

fc=o,+ - 

where the eigenvalues A& are smooth on M 2 \{0} and homogeneous of degree 1, with values 

A (r / ,C) = 0, \ + (r 1 X) = V^f + e, A_(ry, C) = ~^F+C, (63) 

and the projectors pk are smooth on M 2 \{0} and homogeneous of degree 0, with the following values 

Po is the orthogonal projector onto Span I I „ 

p± is the orthogonal projector onto Span 



where Z = I 77 Z~ 



Using Fourier transform, this allows to state the 

Lemma 4.9. Take u m G #°°(R 2 ,C 6 ). 

Then, the unique solution u G C(R, H°°) to the Cauchy problem 

d T u + M 2 (0, d y , d z )u = 0, M| T=0 = u m , 
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is given by 

u= ^ Uk, 

fc=0,+ - 

where each u k = p k (D y ,D z )u is characterized by 

(d T + i\k(D y , V z ))u k = 0, u fc | T=0 = Pk(D y , U z )u m . 

In passing, Lemma 14.91 implies equation (1551) induces for u{] the splitting 

u o _ o , o , o 

u ~~ u 0,0 ' u 0,+ ' u 0,-> 

with 

u° ok =p k (D y ,D z )u° , (d T + i\ k (D y ,D z ))u° Oik = (fc = 0,+,-). (64) 

Lemma [4.91 also implies that equations (1501) . (l51~j) . and (I58|) have the form 

(0r + iAft(D v , D z )) uJJ = • • • (A; = 0, +, -), 
(<9 T + zAo(A,, !>,)) nu^° c = • • • , (<9 T + iX (D y , D z )) Nj'° = • • • , 

where the various right-hand-sides only depend on u° and N°. 

With this observation in mind, the next step consists in analyzing the above equations using 
the average operators introduced in [30]. They allow to describe nonlinear interactions between the 
various modes 0, +, — in the equations at hand. 

For each k — 0, +, — , and u G C([0, £*] 4 xIR T , H°°(Rl ^xT 2 ?)), we define (omitting the dependence 
upon t, x and 9) 

G s k u(T, y, z)--=\J ?~ X {e tsX ^u(T + s, rj, C)) ds, 
with J 7 the Fourier transform in variables y, z. We also define the limit (if it exists), 

G k u(T, y, z) := lim G%u(T : y,z). 

S— >+oo 

The average operator G& performs the average along the bicharacteristic curves of the operator 
(<9t + iX k (D y , D z )). Naturally, G coincides with the usual average with respect to T, due to A = 0. 
The following properties of the average operators G k are useful. 

Proposition 4.10 (borrowed from [30]). Let k G {0,+,—}. 

(%) Ifue C([0,£*] t x R T ,H°°(Rl tyiZ x TTf*)) satisfied {d T + iX k (D y ,D z ))u = 0, then G k u = u. 

(ii) If f e C([0,t*] t x R T ,i?°°(R3^ 2 x T 2d )) satisfies G k f = 0, then any solution u to (d T + 
iX k (D y ,D z ))u — f is T-sublinear. 

(m) Ifu G ^([O^^xRt, H°°(R 3 x ^ z xT 2 e d )) isT-sublmear, then we have G k {d T + iX k (D y , D z ))u = 0. 

(w) Let a collection (u e ) < £ < L C C([0,t*] t x R T ,H°°(R 3 XjytZ x T 2d )) satisfy {d T + iX h (D y , D z ))u e = 
with ki G {0, +, — } for all £, and set u := u - ■ -Ul- 

Then, if ki = k for all £, we have G k u = u, else, if kp^ k for some index £, we have G k u = 0. 



2 Here, solutions to equations of the form (dr + i^k(D y , D z ))u = / are meant in the mild sense, see remark 14.131 
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Our last step consists in applying all above considerations. Indeed, Proposition 14.101 (ii)-(iii) 
asserts that equation (dr + i^k(D y , D 2 ))w = / possesses a T-sublinear solution if and only if G k f = 0, 
while Proposition 14.101 (iv) allows to explicitely compute G k f when / is a product of solutions to 
{d T + i\ ki {D y ,D z ))u e = 0. 

Concerning equation (150]) , with the wave structure given by (1551) . we conclude that equation (1501) 
possesses a T-sublinear solution u ' if and only if 

p k {D y , D g )Mi(dt, d x ) p k {D y , D,)u° * = {k = 0, +, -). 

These three systems in fact reduce to the trivial condition 

d t u° 0fi = 0, d t u° o>+ = 0, 3tu°_ = 0, (65) 

thanks to the following Lemma whose proof is a straightforward computation. 

Lemma 4.11. Take k e {0, +, -} and (rj, () G 1R 2 \{0}. Then, we have 

Pk(vX)A x p k (vX) = o. 

Concerning equation ( 15T1) . the right-hand-side of this equation only involves (products of) solu- 
tions to [8t + iXk t (D v ,'Dz))'u>e — 0, amongst which only the parts u° sc , u[j and N° are associated 
with the characteristic speed A = 0, thanks to relations (139]) and (14T|) . and thanks to the definition 
of Uq in (164")) . Therefore, equation (15T1) possesses a T-sublinear solution u^° if and only if 

(d t + v(D fl ) d x ) u° sc - i a(D fl ) A,, z u° sc = (66) 

n(o,^Tr (r(z fi 7 (* fi 7 -* • deX 1 - 1) [(jg£ + e° osc ) • r,N°] osc ) ). 

This equation may be transformed further. Indeed, when a G Co, we know from Lemma 14.51 and 
Lemma 14.61 that v(a) = and a(a) = 0. Besides, we already know that the populations N° only 
carry temporal oscillations, according to (137|) . As a consequence, we recover by a direct computation 
that equation ( |66l) implies 

V«GC , <9 4 u°=0. (».e. <9 t u° pace = 0). (67) 

(This comes from the fact that the factor i f2 7 (i f2 7 — k ■ 9q x )~ 1 — 1) vanishes when acting on a 
frequency a G Co). When a G C± at variance, equation (1551) provides, using Lemma 1431 and Lemma 
6j the relation 

VaeC ± , (d t ±d x )u° a ^7^—&y,z<= (68) 



n(o,*Tr (r(fi 7 (fi 7 -* • aO" 1 - 1) [(E° ofi + E° osc ) • r,N°] ( 



Lastly, concerning equation ( 1581) . the right-hand-side involves in the similar fashion only products 
of the profiles N°, Tv°, and C 1 ' , which, considering the relation (14"3]) . reduces to products of the 
profiles N° and _E°. Amongst these profiles, only the parts N°, E® sc , and Tvqo are associated with 
the characteristic speed A = 0. Therefore, equation ({58"]) possesses a T-sublinear solution N^° if and 
only if (here we plug relation (143]) ) 

Va G C U {0}, 
d t N° a = W%N° a - 



{e% + e o sc ) • r, (i n, -k ■ d dl y l [(Ko + e° osc ) ■ r, N°] od . (69) 

J d,a 
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Eventually, we have now obtained the set of equations (j55|) . floTl) . (J5S]) . and (J55|) as a set of 
necessary and sufficient conditions to be able to find T-sublinear solutions to ( |50|) . (!5T|) . and (|58|) . 
respectively. This completes our effort in finding a system that completely determines the dominant 
profiles u° and N°. Note in passing that equation (168|) in fact reduces to 

\faeC ± , (d t ±d x )u° a ^i-J—A yjZ u° a = (70) 



u{o,iTr(r(n 1 (n 1 -k-a 1 )- 1 -i)[E^ ime -r,N ] a )). 



due to the fact that N° = whenever a ^ C U {0}. 

In conclusion, we have recovered equations (El, ([26]), ([27]), ([28]), ([29]), ([30]), §2) that had been 
announced in the introductory part of this paper. 

4.3.2 Analysis of coherences: exponential growth 

The above analysis is not possible in the case of f lBT]) : equation (T4"6"]) . ruling the evolution of C ' 1 as 
a function of T, is not a constant coefficient system. Instead, we consider it as a non-autonomous 
system of linear ODE's, parametrized by t, x, y, z, 8, with smooth and bounded coefficients. This 
point of view provides the simple 

Lemma 4.12. Let (u°, N°) e C([0, tJ t x R T , H°°(R 3 X ^ Z x T 2d )) be a mild solution to ([38]). (1391). (|47D. 
dSSJ, flSS}, dSSD- T/ien, 

(%) (u°,N°) G C£°([0, £*] t x R T , i7°°(R| x T 2 e d )) is smooth and bounded uniformly with respect to 
all variables, as well as all its derivatives. 

(ii) associated with these values of (u°,N°), take any solution C 0,1 to (j4"6"l) with an initial data 
belonging to C°° ([0, £*] t , H°°(Rl yz x Tj d )) . Then, for any multi-index /x G N 5+M , there are 

x,y,z 

\d^ y , z ,eC ' l (T)\<K ie ^ T . 



constants Ki,K 2 > suc/i t/iat ; uniformly on [0,t+] t x M T x ]R 3 x T^ d ; we have 



Remark 4.13. All equations that are referred to here are either of the form dtv = f(v,Dv,D 2 v), 
or of the form Otv = f(v,Dv,D 2 v), where v is assumed to have H°° smoothness in (x,y,z,6), 
where the symbol D means differentiation with respect to (x,y,z,6), and f is a possibly non-linear 
function that depends on the equation. The notion of mild solution is then the standard one : we 
mean a solution to the integral equation v(t) = v\t=o + L f(v,Dv,D 2 v)(s) ds or to v(T) = v\t=o + 
fif(v,Dv,D 2 v)(S)dS. 

Proof. Boundedness and smoothness of (u°,N°) with respect to (x,y,z,8) is obvious. Concerning 
the variables t and T, we exploit the fact that t belongs to a compact set. We also exploit the structure 
of the relevant equations. More precisely, u[j satisfies (<9y + M 2 (0, d y , d z )) u[j = and <9 t u[j = 0. This, 
together with Lemma H~9l provides boundedness and smoothness of u° with respect to t and T. Next, 
N° satisfies <9yN° = together with equation (169]) . an equation of the form c^N = /(u°, N°) where 
/ is smooth. This provides boundedness and smoothness of N° with respect to t and T. We have used 
here the Diophantine Hypothesis 13.31 to conclude that / is smooth. Lastly, u° sc satisfies <9*ru° sc = 
together with equation (1661) . an equation of the form <9iU° sc = / (u° sc , D e u® sc , Dju^ sc , u°, N°), where 
/ is smooth. This provides boundedness and smoothness of u° sc with respect to t and T. 

Point (ii) now comes from the Gronwall lemma, together with the fact that equation (j4"6l) has 
the form <9tC 0,1 = / (E°, C ' 1 ) where / is bilinear. □ 
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Remark 4.14. Under the assumptions of the above Lemma, pushing the analysis further, we may 
consider a solution C 1 ' 1 to ( IS"?]) , and look for the available estimates on C 1,1 . This function is 
associated with a solution N 1,1 to (148]) . In view of the above Lemma, the latter obviously satisfies an 
exponential bound of the form \d^ Tx 6/ N 1 ' 1 (T)| < Kie K2T , where K\ and K 2 depend on a but not 
on T. Hence C 1,1 satisfies similarly \d^ Tx e C 1 ' 1 (T)\ < K 1 e K2T . As a consequence, we recover 

\^/eC 1 ' l (t,x,y,z,t/Ve,9)e-^ t/e \ <K x ^~ee {K ^ ll ^ )tl ^ = 0^e). (71) 

4.4 Solving the profile equations 

4.4.1 Computing the dominant profile U°. 

The first profile U° = (u°,C°,N ) is constrained by the polarization conditions (135]) . ( 1361) . ( 1371) . 
Besides, its components are propagated in various ways. The average Uq does not depend on the 
slow time t (equation (165]) ). and it satisfies the linear hyperbolic equation [d T + M 2 (0, d y , d z )) u° = 
with respect to the intermediate time T (equation (164]) ). The oscillating part u° sc , and the populations 
N°, do not depend on T (equations ( 139]) . (147])). and they satisfy nonlinear evolution equations with 
respect to time t (equations ( 166]) . (169])). Coherences C° are only constrained to satisfy the nonlinear 
ODE ( 146]) in time T, and the slow time t only acts as a parameter here. 
We are now in position to state the 

Theorem 4.15. Let s > (3 + 2d)/2. Take a function^ 

H° = U°(x,y,z,9 ), 

which belongs to H S (R X x Tg Q ). Assume that the component C° of U° satisfies the following 
polarization condition 

V(m, n) G 7? d satisfying (m, n, a) <£ lZ(k) for all a G Z 2d , we have C^ n = 0. (72) 

Then, there is t* > and a unique function 

U° = (u^C^e-^N ) e C([0,tJ t x [0,+oo[ T x[0,+oo[ CT ,/F(l^ j2 x Tg, A) )), 

solution to the polarization conditions (135]) . (136]) . (l37j) . and which satisfies equations (164]) . (165]) /or 
the average field u[j ; equations ( 139]) . ( 166]) /or i/ie oscillatory field u° sc , equations ( 147]) . (169]) /or i/ie 
populations N° ; and equation ( 146]) /or t/ie coherences C ' 1 . Uniqueness and existence is guaranteed 
provided we impose U° satisfies besides the two initial constraints 

Ul_ t _ e _ = U°, o^c? C 0,1 | T=0 zs independent oft. (73) 

7/m addition we take an integer £ G N and assume that s > (3+2d) /2+£(a+2) , where the positive 
real a is the one entering the Diophantine Hypothesis \3. 31 then we recover the higher regularity 

U° G C e ([0,a x [0,+oo[ T x[0,+oo[ (J ,ir-^ +2 > (R^ x T^ oA) )) . 



3 The reader should be cautious about the fact that the initial datum only depends on 8q, and not on 9\. Recall 
that eventually ^o will be replaced by kx/e while 6\ takes the value —kt/e. 

4 Due to the fact that the coherences only need to satisfy an evolution equation in the intermediate time T, note 
that the slow time t needs to be treated separately here, including in terms of initial data. Note also that the initial 
data is not decomposed into modes, contrary to the solution itself. 
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Proof. 

First step : identifying the initial values of the profiles. 

The first and most important step consists in understanding the polarization conditions, and how 
they generate the relevant initial data for the different modes. The data U°(:r, y, z, 9q) is the value 
of the solution U 0,K att = T = a = and 9\ = 0. We prove that it defines the value of U°' K at 
t = T = for all values of 6 1 eT d . 

Firstly, concerning the fields, the polarization constraint fl35|) provides, for the (to be defined) 
solution u°, the relations u° = u 0,0 and Ilu = u°. On the other hand, we may write down the 
Fourier transform of u° in the sole variable 6q, together with the Fourier transform of u° in (9q, 9\), 
and obtain 



u°(x,y,z,9 ) = J2^(x,y,z)e^ 6 



Y, u° a (t,T,x,y,z)e ia < e °M. 



u°(t,T,a,x,y,z,9 ,9 1 ) = 

aeCU{0} 

Taking t = T = a = 0, 9\ = in the second equation, and equating Fourier coefficients, then provides 



V/3 G Z d \ {0}, u° = (u° (m + u^_ p) + u°3 i0) ) 



u. 



u 



t=T=0 



(0,0) 



t=T=0 



where u?*^ = tt + u% ^, and u?o m = ■n u? /3 _ /3 x, and u^qn = 7r°u^ - ) , thanks to the polarization 
conditions and using the value of II (see ( |34"|) ). Hence, applying successively the three operators n + , 
7r~ , 7T° on both sides of the first equality, and using the obvious orthogonality relations ti + it~ = 
and so on, we recover the necessary relations 



U (0,0) 



u 



t=T=0 





Ho> 



(P,0) 



t=T=0 



vr+u°, 



and, when (3 ^ 0, 

o 



u 



(P-P) 



t=T=0 



k R% 



u 



08,0) 



t=T=0 



7r°uJ. 



This terminates the analysis of the initial conditions for the fields. 

Secondly, concerning the populations, the polarization constraint (j57|) provides the simpler re- 
lations N° = N 0,0 and N° = whenever a £ Co U {0}. Writing down the Fourier transforms of 
N° and N° as in the previous paragraph then provides N (x, y, z, 9n) = \ No(i, y, z)e ll3 ' d °, and 

N°(t,T, a, x,y, Mo A) = ^ eZd N°(i,T,x,y,2)e^°. Specifying t = T = a = 0, 9 X = 0, and 
equating Fourier coefficients, gives the simple value 



V/3eZ a 



N? 



t=T=0 



m- 



This terminates the analysis of the initial conditions for the populations. 

Lastly, concerning the coherences, the polarization conditions ( ]36l) asserts C^ n a = whenever 
(m,n,a) <£ lZ(k). Writing down the natural expansions of C° = C ' 1 e~ a and C°, we recover for 
each value of m and n the identities C^ n (x,y, z,9 ) = Yl,Rei d — m,n,p{ x iVi z)e l ^' 9 °, together with 
C° m Jt,T,a,x,y,z,9 ,9 1 ) = (£ Q6Z2d C°^ a {t i T,x i y i z)e ia < e °^) e~°.' Specifying t = T = a = 0, 
9\ = 0, gives using the polarization conditions 



V/3 e Z d , C ^ = J2 C° m ) 



n,[p,ai) 



(74) 



ai£Z d 



Now, two cases occur, depending on the value of (m, n) 
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• First case : for any a G Z 2d , we have (m, n, a) (£ TZ(k). 

In that case, equation (174)) is automatically satisfied since the right-hand-side vanishes due to the 
polarization constraint on C 0,1 , while the left-hand-side is assumed to vanish thanks to the additional 
constraint (1721) we have set on the initial function C°. 

• Second case : there exists an a G Z 2d , such that (m,n,a) G lZ{k). 

In that case, the solution ol\ G U 1 to the equation u(m,n) — k ■ ot,\ is unique, thanks to the 
Q-independence of the coordinates of the wave-vector k. Recall that this equation defines the set 
IZ(k). Therefore, given any (3 G Z d , equation (174"|) reduces to 

C ' 1 „,~s = C° B , andVax^o?!, C ' 1 IR > = 0. 

This terminates the analysis of the initial conditions for the coherences. 

Second step : solving the evolution equations. 

We now consider the evolution problem with respect to the times t and T, with the above derived 
initial data. 

The field average Uq decouples and may be determined first starting from the initial value ufji _ _ , 
by using (d T + M 2 (0, d y , d z ))u° = and d t u° = 0. 

Next, we may solve the equations on (Uq SC ,N°). The evolution with respect to T is trivial since 
<9 T (u° sc , N°) = 0. There remains to solve the coupled nonlinear system ( 15S|) . (155|) . an evolution 
equation in t. To do so, we use a standard iterative scheme (see for example [I]), and introduce the 
iteration 

(d t + y(De) d x ) u^ +1 ) - ia(D fl ) A,,,u^ +1 ) = 

n(0,zTr (1X0,(0, -fe-cy- 1 - 1)[E°J^ -r,N°'W] osc ) 



a N 0,(n+1) = ^ u N 0,(n+1) 



(e ^ + < ) • r > (« ^ -* • SO" 1 ! W° + £ °o) • r, N°-W] od 

(3S£° + ^o°o) ■ r, (i O, -A; ■ fy,)- 1 ^"' + £ ° ) ■ r, N°-W] od 



d, space 
Jd,o' 



with initial data f(u° sc ), _ t _ , N°| t=T=0 ), and an initial "guess" set to, say, (u OS c ,N 0, (°)) = const = 

(( u osc)l_ t _ , N°| t=r=a ) for any t (and T). For this linearized scheme, usual energy estimates are 
available in any Sobolev space H s , to which the skew-symmetric operators v(D#) d x and i&(Dg)A ytZ 
do not contribute. They are 

po,(m+i)^ _ ^.(n+ij^n < /" e C(t-,) || F (C7<',M)( 5 ) _ F(U°' {n) )(s)\\ H3 ds, 

Jo 

where F{U ) stands for the right-hand-side in the above iteration. Now, standard nonlinear tools 
assert that the function F acts in a locally Lipschitz fashion on H s = H S (M. 3 x T 2d ) provided 
s > (3 + 2d)/2. The existence and uniqueness of a solution (Uq SC ,N°) to ( |66|) . (169)) on some time 
interval [0, t*\ then follows. 

Concerning the coherences C 0,1 , the statement of our Theorem imposes to choose them indepen- 
dent of t at time T = 0. This constraint is only a (pratical) way to fix the value of C 0,1 on the set 
{T = 0} (recall that C ' 1 satisfies an evolution equation in time T only). Equation ( H6|) then asserts 

d T C^ 1 = t[E°-T,C ' 1 ]. 
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This equation is enough to uniquely determine C ' 1 for any t and T, starting from its known values 
on {T = 0}. 

The polarization conditions (|35|) . (|36|) . (|37|) commute with equations (164"]) . (165]) for the average 
field u[], with equations (|39|) . (|66|) for the oscillatory field u° sc , with equations ( 147"]) . ( 169]) for the 
populations N°, and with equation (T4"6"]) for the coherences C ' 1 . Hence by uniqueness of the solutions 
to these propagation equations, polarizations are preserved along the evolution. 

There remains to study the higher regularity of the solutions we have exhibited. The equations 
on the mean field, the equations on the populations, and the equations on the coherences clearly 
do not induce any loss of smoothness, i.e. provided u° sc is H°°, the derivatives <9 t (u[], C ' 1 , N°) 
and <9t(uq, C 0,1 , N°) have the same H s smoothness as (u(j, C 0,1 , N°). The loss of smoothness comes 
from the Schrodinger-like equation ( 16"B"]) on u° sc . When differentiating this equation with respect 

to t indeed, an additional factor a,(a)/S. vz = ±— ; A v , comes into play which causes, from the 

Diophantine Hypothesis 13.3] a loss of Sq = a + 2 derivatives (two derivatives are lost due to the 
Laplacian, while a derivatives are lost due to small divisors, seen in Fourier space). This argument 
explains why U° G C 1 (if s_s °) and, more generally, why U° G C £ (H s ~ es °). □ 

4.4.2 Computing the correctors U 1 and U 2 . 

The first corrector U 1 = (u 1 , C 1 , N 1 ) is built up in the following way. On the one hand, we have 

u 1 = u 1 ' + u 1 ' 1 e~% C 1 = C 1 ' + C 1 ' 1 e-% N 1 = (Nj'° + N^J + N 1 ' 1 e" CT , (75) 

thanks to (PD , {HP , PJ . Besides, the contributions (1 - II)u£°, u 1 - 1 , C 1 - , C# n>a l( m , n ,a) W ), 
N 1,1 , are prescribed as functions of the dominant profile U° through the polarization conditions ( 1401) . 
( ]4il . l j4"3"|) . (jloTl . (j4"8~l) . respectively. They are thus known and smooth functions of U°. Let us write 
these relations symbolically as 

((1 - HK'° C , u 1 ' 1 , C 1 ' , C^ l (m ,n,a )m k), N 1 - 1 ) = (U°) . (76) 

The remaining terms 

i-i 1 >° nnM n 1 - 1 i ivt 1 ' _i_ Tvi.o 

U > iiU osc' ^m,n,a 1 (m,n,a)e^(fc), ^o + iN space' 

are then determined through linear evolution equations in T, namely through ( 1301) . ( T3TT) . ( 1571) . ( |58l) . 
respectively. These equations read, after some simplifications, 

1,0 



(3r + M 2 (0, d y , d z )) < u = -A^itf , «9 T nu^ c = 



irz^o 



VogCqU {0}, d T Wf = - E u -T,(i fi 7 -k ■ dg )- l [E" ■ T, N 



od 



+ 

V(m,n,a)eH(k), d T C l m \ a = -d t C Q + + i [E°'° ■ T, C 1 ' 1 N J# — , 



d,« 



d.Q: 

(£ o sc + E° 0fi ) • T, (i fi 7 -k ■ de^KE^ + E° 0fi ) ■ T, N°] o ; 
_Ano,i + i [ E o,o . p ; c i,i + N i,i] 

+ i [E lfi • T, C ' 1 ]^ + i [E 1 ' 1 ■ T, N°'°] m , n , Q . (77) 



Note that the right-hand-side of the above system only involves known functions of U° at this stage. 
In order to solve this system, there remains to impose as we did for C ' 1 , 

(uo' ,nu^,N^ l QeCoU{ o},C^ nQ l (mi „ ia ) e ^ (fc )) _ is independent of t. (78) 
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This completely determines U 1 as a function of U°. 

The second corrector U 2 = (u 2 , C 2 , N 2 ) is built up in the following way. On the one hand, we 
have 

u 2 = u 2 <° + u 2 ' 1 e'\ C 2 = C 2 ' + C 2 ' 1 e~° + C 2 ' 2 e" 2CT , N 2 = N 2 ' + N 2 ' 1 e~° + N 2 < 2 e" 2CT , (79) 

thanks to (IMl) . (joT)]) . (|60l). Besides, the contributions (l-n)u 2 £, u 2 ' 1 , C 2 ' , C% n>a l( m , n ,a)£R.(k), C 2 ' 2 , 
N 2,0 l a ^c u{o}? N 2 ' 1 , N 2,2 , are prescribed as functions on the two first profiles U° and U 1 through 
the polarization conditions (JS2D? dSS]), (l5"5"l) . ( ISljj) . (1B3|) . (jSHJ), (EOD, dSH]), respectively. They are thus 
known, smooth functions of U° and U 1 . Some parts of the corrector are free and may be chosen 
arbitrarily, namely u 2 ,' , Ilu 2 £, C 2 ^^ l( m ,n, a )eit(k), N 2,0 l Q6 c u{o}- We make the most simple choice 
to set these contributions to zero. Let us write all these relations symbolically as 



u 2 = ^(u°,u 1 ). 



(80) 



All these considerations, in conjunction with Theorem 14.151 lead to the 

Proposition 4.16. For all leN, there exists s(£,a) such that the following holds. 

Take an s > s + (3 + 2d)/2. Take an initial data U° in H S (M^, x Ti) which satisfies the 



constraint (IT2j) . Take the profile U° associated with these data through Theorem J^.15\ Take two 
initial data U 1 in H s (Rl >y>z x Tg o ) and U 2 in H s (Rl,.„ x 



,. ,, . ,, Tg () ). Then, there are unique correctors 



u i,o + u i,i e -<r 



U 1 - i C^ + C 1 '^ 



U 2 



u 2,0 + U 2,l e - CT 

C 2,0 + C 2,l e -<x + C 2,2 e -2 CT 

N 2,0 + N 2,l e -<x + N 2,2 e -2o 



in C e ([0, t*] t x [0, +oo[ r x [0, +oo[ CT , H* 



-s{to3 



x,y,z 



T 2 g 9l ))), which satisfy the constraints (ITS]) . (1781) fi 



or 



the first corrector, the constraints flHDj) /or t/ie second corrector, and the evolution equation (158]) for 
the first corrector, together with the initial constraint 



u 1 


T=t= 


= u 1 u 2 

:o-=0,6»i=0 


T=t=o-=0,6»i = 


= U 2 . 

=0 


Proof. The initial constraint U 1 


= U 1 reads 

T=t=cr=0,6»i=0 


(u^ + u 1 ' 1 , 


c i,c 


+ C^,(Nr + N s 1 p ° c 


J+N 1 ' 1 ) 


T=t=O,0i = 



u 1 . 



This, together with (TTB|) written at time T = t = 0, prescribes the value of 



((1 — II)u ^ c ,u ' ,C ' , C n ' jm>Q l( mj „ ia )^(fe),N ' ) 
Hence, taking the difference, we recover the value of 



T=t=0 



(U°) . 



U + iiU osc' 2^^ ■n,m,a 1 (rn,n,a)GTl(k)e ,JN + iN space 



T=i=O,0i=O 
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an explicit, linear function of U° and U 1 . Fourier transforming U° and U 1 in the variable 9 , and 
Fourier transforming 



U Q -t- llU osc , / j ^ n ,m,a - L (m,n,q)g^(fc) C , V\q T 1"N 



1,0 

space 



in the variable (#o;$i) then allows to deduce, as we did in the proof of Theorem 14. 15[ the value of 
(uj'°, n<°, C^ na „ N$) I _ _ , whenever a G Z M \ {0}, (m, n, a') G K(k), and a" G C U {0}. The 

requirement (1751) then gives the value of (u ' , IIu^ , C^ n a ,, N a „) on the whole set {T = 0}U{£ = 0}. 
Equation (J58j) in turn provides the value of (u ' , IIu^ , C^ n Q ,, N tt '„) for any value of T and t. The 
remaining part of U 1 , namely the value of ((1 - II)u^°, u 1 - 1 , C 1 ' , C^ ma l( m , n , a )$ii(k), N 1 ' 1 ) for all 
values of T and £, is next given by (f> (U°). This terminates the computation of U 1 for all values of 
T and t. The relation U 2 = ip (U°, U 1 ), see (15U1) . prescribes U 2 for al values of T and t. 

Let us come to regularity issues. Thanks to the constraint (176"]) . the C £ (iJ s-s (^ a )) smoothness 
of (1 - IIK£, u 1 ' 1 , C 1 - , CJi^l^a)^), N 1 ' 1 is simply the consequence of the C e (H s -' s ^) 
smoothness of U°. For Uq' , Ehi^g, N^° l QeCo u{o} and C^ Q l(m,»,a)ew(k)» the C°(F s ~ s(0 ' a) ) regularity 
stems from the classical smoothness of solutions to linear hyperbolic systems with source terms lying 
in C°(H s ~ s (°' a ^). Note that in the last equation in (j58|) . the regularity of the source term <9 t C 0,1 is 
a simple consequence of the regularity of solutions to ODE's depending on a parameter (here t). 
Note also that the overall loss 3(0, a) comes from both the involved derivatives of U° appearing in 
the source terms, and from small divisors (u(m,n) — k ■ ai)^ 1 and M\(—k ■ «i, k ■ «o) _1 acting on 
U°, that appear in the source term as well. Differentiating (177)) and (175)) with respect to T and t, 
and applying the same argument, eventually provides the C t (H s ~ s ^ ,a ^) smoothness of U 1 . Relation 
U 2 = tp (U°, U 1 ) provides the C e (H s -~ s ^) smoothness of U 2 . □ 

Remark 4.17. One can prove that for any £ we have U 1 G C £ (_£p- max ( 3a + 1 > a + 2 M( a + 2 )) ; and U 2 G 

nil rrs-max(4a+2,3a+3)-^(a+2) \ 

5 Convergence 

5.1 The residual 

Writing down the full profile representation of the residual 



l -d -1 



r £ (t, x,y,z) = L[ d t , d x , —d y , -=d z U £ app - F £ (U app ) 



we have, using Proposition | 

r £ (t, x, y, z) = TZ £ {t } x, y, z, T, a, 



T=t/ y /e,<T='yt/e,(eo,ei)=(kx,-kt)/e 



3 

with TZ £ = Y^ YVe 1 r j ' K {t,x,y,z,T,e)exp{-Ka). 

j=-2 k>0 
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We now find how r 1 , r 2 and r 3 depend on the profiles U 1 , U 2 and U 3 . Symbolically, introducing 
first order differential operators L\ and L 2 , a matrix Lq, and some bounded, symmetric, bilinear 
operators Be, -B(c,n) (with values on the E- and (C, N)-components, respectively), we have 

r 1 = L 2 (d T , d y , 8 Z )\J 2 + L x (d t , d x )\J l + 2B E (U°, U 1 ) + 2B (C>N) (U°, U 2 ) + ^(^(U 1 , U 1 ) + L U\ 
r 2 = Li(^, ^)U 2 + 2B E (U°, U 2 ) + £ B (U\ U 1 ) + 2£? (c ,n)(U 1 , U 2 ) + L U 2 , 
r 3 = 2I^(U 1 ,U 2 ) + £ ( c,N)(U 2 ,U 2 ). 



Naturally, once C°° profiles have been built thanks to Theorem 14.151 and Proposition 14.161 regularity 
of these residual profiles is clear. The question we now need to face is to evaluate the size of these 
residuals as the intermediate time T grows unboundedly. 

Concerning the dominant profile U°, the components u° and N° are bounded uniformly in T, while 
the component C° has exponential growth, i.e. C° and all its derivatives have size K\ exp(K2T) as 
T increases, see Lemma [4.121 

Concerning the first corrector U 1 , the components u ' , IIu^ and N ' + Ng'° ce are precisely 
constructed so as to be sublinear in T, see Section 14.31 This is the key point. On top of that, 
the component (1 — n)u^° is bounded, thanks to ( l4"Uj) . the component u 1 ' 1 has exponential growth, 
thanks to (141 p . the component N 1,1 has exponential growth, thanks to (1481) . the component C 1,0 is 
bounded, thanks to ([43]) . and the component C 1 ' 1 has exponential growth, thanks to (|4"5"|) and (1571) 
in conjunction with the Gronwall Lemma. 

Concerning the second corrector U 2 , we know from (I8UI) that U 2 is an explicit, linear function 
of U° and U 1 , functions that have at most exponential growth. Hence U 2 has at most exponential 
growth in T. The only difficulty may then come from the component U 2 ' , whose exponential growth 
will not be eventually compensated by a decaying term e~ a or so. In that direction, we observe that 
u ' , LTu 2 ^, and N 2,0 l Qe c u{o} are conventionally chosen to vanish. On the other hand, the component 
(1 — Il)u 2 ^ is bounded thanks to (15"2l and to the boundedness of (u°, N°, u^), the component C 2 ' 
is bounded thanks to (15T)1) . to the boundedness of C 1 ' , u°, N°, and to the sublinearity of N 1,0 , while 
the component N 2,0 l a <£c Q u{o} is bounded thanks to (1511 . to the boundedness of C 1,0 , u°, N°, and to 
the sublinearity of N 1 ' . 

As a conclusion, we have now established the 



Lemma 5.1. Given the C°° profiles provided by Theorem \4.15\ and Proposition \4-lo] with the choice 
s = +oo ; the following result holds. 

For allfi G N 5 + M ; there are constants K 1 , K 2 > such that, uniformly on [0, £*]£xR T xIR 3 ^xT 2 /^ 
for all k G N* ; j = 1, 2, 3 ; we have 



\ u t,T,x,y,z 



r j ' K (T)\ <K ie K2T , and ^|a M TT „ ? ,r J '°(T)| — > 0. 



As a consequence, we get for the residual 



V/x e N 3+d , 



sup 

te[o,t*; 



^x,y,z,e ^ £ \T=t/y/E,<T=-yt/efir- 



-kt/e 



L 2 e- 

x,y ,z 



5.2 Stability 

In this section, we prove our main Theorem. 
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Theorem 5.2. Let the profiles U° ; U 1 , U 2 G C°°([0, £*] x [0, +oo[ 2 , #°°(M 3 x T 2d )) be given by Theo- 
rem \4-15\ and Proposition \4-lb] They provide us with the approximate solution U app = U| pp (£, x, y, z) 



given by equation (|33|) . 

Then, for any s > (3 + 2d)/2 and any familly (5 e ) £>0 C H S (R 3 x T d ) such that \\5 £ \\h s — ^0 ; there 

e— >0 

is Sq > such that for e g]0, £o], the Cauchy problem 

l (%&, 4=9 V » 4=^ ue = F£ ( ue )' 

Uf t=0 = U| pp|t=Q + 5 £ (x, y, 2, fcx/e), 

/or the Maxwell-Bloch system fl9]) /ias a unique (mild) solution XJ e which belongs to C([0, £*], H S (M. 3 )). 
Besides, for all /i6N 3 such that s — |/x| > (3 + 2d)/2, we have 

\\ d x, y ,z (U £ - U app) IU»([o,t*]xR») ^°- 

Proof. Standard results for symmetric hyperbolic systems ensure that, for e > fixed, a unique 
mild solution U £ G C([0, t e ], i7 s (M 3 )) exists for some t e > 0. The difficulty lies in bounding t £ from 
below. 

We use a singular system method (cf. [25J), and look for U £ under the form of a profile, namely 
look for U e such that 

U e (t, x, y, z) = U £ (t, x, y, z, kx/e), 

where U e = U £ (t, x, y, z, 6q) corresponds to the following initial data, which is non-singular in e (this 
is a key point) 

W £ (0, x, y, z, )=(j2 Ve~ JlJj (t, x, y, z, T, a, 9 , 0i) + 6%t, x, y, z, 9 ) 

\j=o / I 

\t-T=cr-0,e 1 =0 

It is then sufficient, for U £ to be a solution to (JHJ), that W satisfies 



L (d t , d x + ^d 6o , -^=d y , -^ U e = F £ (W). 



To go on with the analysis, we now set 



Wfpp(t,a;,j/,3,0o) = ^2V^ u3 (^x,y,z,T,a,9)\ Ts=t 



=t/ v / E,<T=7t/e,e 1 =*!t/e 

j=0 



Note in passing that the at-most-exponential growth in T of the various profiles \J^ K implies that 
the family (U^ ) £ is bounded in C([0, £*], H S (R 3 x T d )). We evaluate the difference 



A £ = W - W &vv . 



In the next few lines, we may sometimes write A e = (A £ , A e c , Af^), refering to the u, C and N 
components of A e , respectively. In any circumstance, we have A e G C([0, t £ ], H S (M. 3 x T d )), its initial 



value is A £ | t=0 = S £ , and A e satisfies 



L (dt, d x + \d eo) -Ld y , ±d M } A £ = F^ pp + A £ ) - F £ (W a £ pp ) - ^ t/ ___ _. (82) 



=t/ x /£,CT=7t/e,e 1 =fet/e 
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Hence standard H s estimates provide 

1 d 
2dt 



--||A £ ||^ < (F^W app + A £ ) - F^{UtJ\^) HS + \\n*\\ Hs W\\ Ha . 



Now, the evaluation of the scalar product (F £ (U^ pp + A £ ) — F e (U £ pp )\A e ) H s involves various terms 
that may be ordered in powers of 1/y/s as seen by inspection of the bilinear function F £ in ^. The 
term carrying the weight 1/e is 

- £ {-i [Q 7 A £ c ] - 7 A £ c | A £ ) H s = -1(A £ C | A%. = -1 || A £ c || ^ , 

where the first equality comes from the fact that the operator [Q, .] is skew-symmetric. To estimate 
the other terms, it is useful to keep in mind that U £ is uniformly bounded in H s on the interval [0, £*], 
so there is a constant C, independent of / and e, such that for any t and e we have ||Wf pp (t)|| < C. 
Using this piece of information, the term carrying the weight 1 in (F e (U £ pp + A £ ) — F e (W £ pp )|A £ ) H a 
is clearly bounded by 

n iia £ II 2 

° 11^ \\H S ' 

for some constant C independent of t and e, while the term carrying the weight yfe in is clearly 
bounded by 

v^||A £ ||^, 

The more difficult term is the one carrying the weight 1/v 7 ^- Carefully treating apart all occurences 
of the two terms A e c and C e , all terms being majorized in the most simple fashion, we recover that 
this contribution is upper-bounded by 

C / 

I ||A e II 2 -I- IIA £ II IIA £ II -i- IIP £ II IIA £ II 2 -J- IIP £ I IIA £ II IIA £| 



r- \\\-^C\\H a ~ II "C Hip II" \\H S ~ ll^appllj^s II" \\H S ~ || v '»pp||H s N Cliff II" \\H S 

All in all we have eventually proved, gathering some terms for convenience, 

_||A £ II 2 <--IIA £ ll 2 + — (IIA £ II 2 -)- IIA £ II IIA £ II "l 
jJI^ I Iff — _ ll^cllff + rz Ul cliff + ll^cllfp 11^ \\h b ) 



C 

j llpe II (\\f\ £ || IIA £ II -I- IIA £ II 2 ^ 

/e " a PPli# s Ul^cllff 11^ \\h° + 11^ HffJ 

_l ^ IIA £ II 2 _i_ /? IIA £ II 2 _l \\1? e \\ IIA £ II 

+ O || ZA ||^ 3 + V^ II A || H s ' II ,<J II ff II ^ II ff • 

Hence, there is an e such that for any e < eo we have, 

j -i /~< ||C £ 

IIA £ II 2 < IIA £ II 2 _i_ IIA £ II IIA £ II I (~i II a PPllff ||A £ II 2 -uH-7? £ ll IIA £ II fR^ 

-^||A || H , < -— ||A c || ffs + — ||A C ||^ S ||A \\ HS +G j=. ||A \\ H „ ■-, || A: \\ Hs ||A || ffs . (8dJ 



Now, the two crucial ingredients are 



thanks to Lemma 15.11 , and 



sup \\TZ e \\ Ha -> 0, (84) 



te[o,t ' - l] 



I Cf — i / ill//- 



app 



dt < C y/e, (85) 
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for some C independent of e. This crucial piece of information comes from the fact that 
C e app (t) = C^l^^e"^ + V~eC^\ T=t/v - e e-^ + 0(e), 

which, in conjunction with the at-most-exponential growth of C 0,1 and of C 1 ' 1 (see e.g. Lemma \A. 121 
and estimate (1711) ). provides 



\Ct pp (t)\\ Hs <exp[-—)+Ce, 



t 
Ce 



for some C independent of £ € [0, £*] and e. 

At this stage, an easy argument using the Gronwall Lemma allows to deduce from (1831) . (jHi 
and (1851) . that ||A e ||# 3 is bounded independently of e over the whole interval [0, £*]• Repeting the 
Gronwall argument next shows that ||A e ||# s actually satisfies 

||A e || H . < W\\h» explc f (\LjEpMl + \\K £ \\ HS J dt 7 ) ^>0, uniformly on [0,£,]. (86) 

We skip the whole Gronwall-like argument. Estimate ( 186]) now induces, by Sobolev's injection, L°° 
convergence to zero for profiles A e = U £ — U £ & , and thus L°° convergence to zero for the original 
functions U — V £ , since II U — U*||i«> = \\U £ — W m \\L°° „ . □ 

6 The Transverse Magnetic case 

6.1 The system 

In this section, we present the previous WKB in the particular Transverse Magnetic case, when fields 
take the form 

B= IbA =B(t,x,y), E = I J =E(t,x,y), 

with the additional common assumption (cf. [M]) that the polarization operator T has entries parallel 
to E, namely 

( ° 

Vm,nG{l,...,JV}, r(m,n)= 

\T(m,n) 



Maxwell-Bloch system then reads 

1 

v^ 



d t B e x + -^d y E £ = 0, 



d t B £ y + -^d x E £ = 0, 

y V? 

8 t E £ - d x B £ y + d y B £ x = -^Tr (r fi 7 C e ) - *£ e Tr (T[T, C £ + N £ ]) - v^Tr (r^ Jj N £ ) , 

d t C £ = -- fi 7 C £ + -^E £ [T, C £ + N £ ] od , (87) 

d t N £ = ^=E £ [T,C £ } d + W$N £ . 
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6.2 The Ansatz 



As stressed in Remark 13. 2\ the introduction of an intermediate time T = t/y/e is not necessary here 
(see also Remark 16. II below) . In order to simplify computations, we also restrict here our attention 
to the case of prepared data, which corresponds to the case when 



n =0 =o. 



This second simplification allows us not to use the variable a = yt/e in the sequel. In a nutshell, we 
here consider the simplified Ansatz 



Vl pp (t,x,y) = X) v^U J '(t,x,y,(9)| fl=(fa/ei _ w/e , U j (t,x } y,9) = £ ^Vi(t,x,y) e l 
The characteristic sets Co, C+, C_ and the resonant set TZ(k) are the same as before. 



6.3 WKB expansions 

In this setting, vanishing of the terms r~ 2 , 



r 



1 and r° from Proposition 14.21 reduces to 



' -k-d 9l B° x 
-k-d 01 B° v 



0. 



k ■ d eo E° 
-k ■ d 6l E° - k ■ d do B° y 
(iO y -k-d 0l )C = 0, 

. -k-d 9l N° = 0, 

-k ■ d 9l Bl + d y E° = 0, 
-k-da.Bl-k-d^E 



'ii/ 



0. 



0. 



-k ■ d Bl E l - k ■ de.B] + d y B x 
(i Q, -k ■ deJC 1 = i E°[T, C° + N°] od , 
{ -k-de^^iE^T^C }^ 



(89) 



(90) 



f -k ■ d 9l B 2 x + d y E l + d t B° x = 0, 

-k ■ d 8l B 2 - k ■ d eo E 2 + d t B° y - d x E° = 0, 

-k -d ei E 2 -k- d 9o B 2 y + d y Bl + d t E° - d x B° y =iTi(T fi 7 C 1 ) - i £°Tr (r[T, C° + N ]), (91) 

(z fi 7 -A; • d 6l )C 2 + d t C°'« = % E°[V, C 1 + N 1 ^ + % E X [V, C° + N°] od , 

-k ■ d 8l N 2 + d t N° = i E°[T, C 1 ^ + i E 1 [r, C°] d + W ft N°. 

According to (189]) . the polarization conditions from equations fl35|) - fl37|) become, with the nota- 
tions from Definition 14.71 



£\time 



0. 



Va i C + U C_ U {0}, B° = £° = 



VaeC.U C_ 



5 



,y,tt 



TEl 



C° = 0, 

1 ^ time 



0. 



(92) 

(93) 
(94) 
(95) 
(96) 
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Here, the projector II from Definition 14.31 is given explicitly, and the "prepared data" condition is 
coherent with the vanishing of C°. 

From (I90j) . we get for the average u[] a version of the evolution equation at intermediate scale 
(|38|) under the form 

dyE% = d y B° £fi = 0, 

which leads to polarization consitions, instead of evolution equations, namely 

K = Ko = 0- (97) 

The oscillating part of f]9"U]) produces the following transcription of the polarization condition (|4"U1) 
for the first corrector u 1 

Ma $ C U {0}, u* = 0, 

\/aeC , El = 0, Bl a = -j^dyB^ (98) 

\/aeC ± , Bl a = -^d y El Bl a = T El 

For coherences and populations, equation (143]) is unchanged, 

C 1 = i(in^-k-d 9i y 1 E [T,N°} d , (99) 

and equation (14"8"|) becomes 

NL e = 0. (100) 

Next, the average of fields equations in ([9~T]) is equivalent to 

d y E*=0, i.e. ^0=° ( 101 ) 

(a polarization condition, again, instead of evolution as in equation ( |50l ). 

8 t B% = 0, (102) 

(playing the role of equation ( l65i) ). and 

d y Bl - d x B% = -^(E Tr (r[r, N°])) , 

which reduces to 

d y B l xfi - d x B% = 0, (103) 

thanks to polarizations conditions and spectral properties (193]) . (196]) and (197]) . 

Remark 6.1. Since B° Q does not depend on time t (according to equation (\102\i ), we may impose (as 
in equation (11031) ) that d x By be the y- derivative of an H s function for all times, simply by requiring 
this condition be satisfied at t = 0. But in the general three-dimensional framework, d t By Q = 
is not given a priori, and we need the addition of the intermediate variable T to perform the 
analysis of Section \4-3\ Omitting this intermediate time leads, in the three-dimensional case, to the 
overdetermined (and ill-posed) system ([3"8"]) . fjoU]) 

M 2 (0, d y , d z ) u° = 0, M 2 (0, d y , d z ) ul = -M x {d tl d x ) u°. 
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Oscillations in f[9"Tj) are analyzed as follows. The polarization (j52"|) for u^ sc splits into 

<tnne=(^^)-VL. (104) 

and 



/ 



Va6C ±J 2 z (Sj >a ± £ 2 ) = _L_02£° _ 2 Tr (r ft 7 C£) + z(£°Tr (r ft 7 N°) Q , 

Va^C + UC_U { 0>, ( i |) = ( ,. an)2 : (<; . ai)2 x (105) 

( h- ol\ —k ■ a \ ( 

X V-fc ■ a fc • ai J \d y B l xa + z'Tr (r fi 7 (£°N° - C 1 ) 

where B\ a vanishes for a ^ C (because of ( 198]) ) , and may be chosen arbitrarily when a G Cq . The 
coherence C 1 is given by equation (199]) . Then, the evolution equation (166]) with respect to the slow 
time £ corresponds to 

dt<s P ace = 0, together with d t B% = 0, (106) 

and, using (199]) . to 

1 
ik ■ a\ 



Va e C±, 2(<9 t ± &)££ + - 7I -— — SgJSS = z Tr (r ft 7 C£) - z (£°Tr (r[r, N°]))< 



or (with C* from (199]) ). in other words 

2{d t + y{D e )d x )E° + d^d 2 y E° = % Tr (T(i fi 7 (z fi 7 -A; ■ ^J" 1 - l)(E°[T, N ])), (107) 

(in place of equation ( 151]) ). with v(Dg) given in Lemma [4.51 
For coherences, we get 

C 2 = % (i fi 7 -A; ■ deJ-^E 1 ^, N°] od + £°[I\ C 1 + N 1 ]^), (108) 

and for populations, we have 

N t 2 ime = -»(*• d ei )- l (E°[T, C l ] od ) time: (109) 

and 

a t N° = W t N° + z (£°[r, C 1 ] d ) spacc + z (£°[I\ c 1 ]^, 

i.e. 

d t N° = W i N°- (E°[T, (z fi 7 -A; • ^)- 1 (S [T, N°] od )] d ) space (110) 

- (E°[T, (i fi 7 -k ■ de^iE^T, N°]od)]d)o. 

6.4 Conclusion in the TM case 

The above computations provide us with a set of profile equations leading to a local in time, smooth 
approximate solution U| pp as in 
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• The leading profile (u°,C°,N ) is given by polarizations (J22J, (JM}, flM}, (ESD (C° = 0), (1961), 
as well as (|97|) (for average of fields, instead of an evolution at the intermediate scale). They also 
satisfy an evolution with respect to time t, given by the trivial equations (I102p and (I106p . namely 

®tB yfi = d t B xsp = 0, 

by a nonlinear Schrodinger equation (11071) . and by the Boltzmann equation (IllOp . As quoted in 
Remark 16.11 this a priori overdetermined set of equations is in fact well-posed because of its scalar 
structure. We stress the fact that no rectification occurs for fields at leading order: the only non- 
vanishing average By is constant in time. 

• Correctors are partially determined by polarizations, namely (|98|) . (llOip . (11030 for u 1 , (I9~9l) for C 1 , 
(ITUU1) for N 1 , (TToijl . (ITO for u 2 , (HDHD for C 2 , and (dUD for N 2 . Parts of the correctors that are not 
submitted to these constraints may be chosen equal to zero. 
The Ansatz 

Uapp(£, x, y) = Wf pp (t, x, y, 0)\ e = (k/e ,_ kt/e) 

is consistant with system (IHTj) . We easily prove the 



Lemma 6.2. Given the C°° profiles above, define the residual 1Z £ as in Proposition ^ -4 Then, for 
all fi G N 2+d , there is K > such that: 



sup 

te[o,t*; 



dx )V fiJt £ W=-ktle 



T 2 



<K^l. 



Note that for prepared data, no initial layer is created at leading order, so that we get a 0{^/e) 
estimate instead of the o(l) in Lemma [5.11 

By the same technique as in Section 15.21 we get finally 

Theorem 6.3. Given the smooth profiles above on [0,£*] ; for s > (2 + d)/2 and any familly (5 £ ) £> o C 
H S (R 2 x T d ) such that \\d~ e \\ H s = 0(y/e), there is e > such that for e e]0, e ], the Cauchy problem 
for Transverse Magnetic Maxwell-Bloch system flH7|) . with initial data U| pp| + 5 e (x,y, kx/e), has 

a unique solution U e G C([0, £*], H S (M. 2 )), and for all fi G N 2 such that s — \f/,\ > (2 + d)/2, there is 
K > such that 

\\ d x,y ( U£ ~ U app) IU°°([0,t»]xIR2) < K^/s. 
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